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Abstract 

We calculate the mass spectra and string tensions of SU(2), SU(3), SU(4) and SU(5) gauge 
theories in 2+1 dimensions. We do so by simulating the corresponding lattice theories and 
then extrapolating dimensionless mass ratios to the continuum limit. We find that such mass 
ratios are, to a first approximation, independent of the number of colours, N c , and that the 
remaining dependence can be accurately reproduced by a simple 0(1/N%) correction. This 
provides us with a prediction of these mass ratios for all SU(iV c ) theories in 2+1 dimensions 
and demonstrates that these theories are already 'close' to N c = oo for N c > 2. We find that 
the theory retains a non-zero confining string tension as N c — > oo and that the dimensionful 
coupling g 2 is proportional to 1/N C at large N c , when expressed in units of the dynamical 
length scale of the theory. During the course of these calculations we study in detail the effects 
of including over-relaxation in the Monte Carlo, of using a mean-field improved coupling to 
extrapolate to the continuum limit, and the use of space-time asymmetric lattice actions to 
resolve heavy glueball correlators. 



1 Introduction 



The non-perturbative physics of QCD continues to be largely impervious to analytic attack. 
Thus 't Hooft's proposal to consider SU(iV c ) gauge theories (with quarks) as perturbations in 
powers of 1/N C around N c = oo JIJ. remains of great interest. In many ways the N c = oo 
theory is much simpler than the physically interesting N c = 3 theory, and the fact that the 
phenomenology of the SU(oo) quark-gluon theory appears to be strikingly similar to that 
of (the non-baryonic sector of) QCD Jl|, ||. motivates the suggestion that the physically 
interesting SU(3) theory might be largely understood if we could solve the much simpler 
SU(oo) theory. Unfortunately an analytic solution of the latter still eludes us, even if much 
progress has been made in understanding aspects of its structure f|, ^ ||. 

This situation has motivated a number of computational explorations [0. Almost all of 
these have used the fact that the lattice SU(oo) theory can be re-expressed as a single plaquette 
theory 0. Although these calculations have produced interesting results, the approach suffers 
from a basic problem: it tells us nothing about the corrections to the N c — oo limit and so 
cannot address the critical question of how close SU(3) is, in fact, to SU(oo). The twisted 
Eguchi-Kawai approach also suffers from the fact that the space-time volume described by the 
theory is finite and related to N c . 

In this paper we take a more direct approach to the problem. We calculate the continuum 
properties of SU(2), SU(3), SU(4), ... theories by simulating the corresponding lattice theories. 
We then compare these properties and see how well they can be described by simple corrections 
to a common N c — oo limit. This approach has the advantage that it will tell us just how 
close the physically interesting SU(3) theory is to the simpler SU(oo) theory. It has, of course, 
a potential disadvantage: if the theories we consider are far from the N c = oo limit then we 
will have learned nothing about the physics of that limit. Fortunately, as we shall see, this 
turns out not to be the case. 

Ultimately we would like to consider SU(iV c ) gauge theories coupled to light quarks (in the 
fundamental representation) in 3+1 dimensions. But this is beyond our current computational 
resources. So our first restriction is to disregard the quarks and focus instead on the pure 
gauge theory. Since the non-perturbative physics of QCD is largely driven by the self-coupling 
of the gauge fields this is certainly a physically relevent problem. Moreover one expects the 
pure gauge theory to have leading corrections that are 0(1 /N^) rather than the 0(1/N C ) that 
one expects with quarks. Thus the onset of large- N c physics should be easier to spot. 

Our second restriction is to consider the pure gauge theory in 2+1 rather than in 3+1 
dimensions. Although it is less obvious that this leaves us with a physically relevant problem, 
we shall argue in the next Section that gauge theories in 3 and 4 dimensions are sufficiently 
similar that this is likely to be so. Moreover it turns out that in D = 2 + 1 we can calculate 
the continuum properties of gauge theories with such accuracy that there is little ambiguity 
in our final conclusions. This is not the case in D=3+l where the preliminary calculations of 
this kind || are very inaccurate in comparison. 

Motivating our calculations are several questions of particular interest. What we know 
about the large-A r c limit of gauge theories essentially comes from considering Feynman dia- 
grams to all orders. Such considerations indicate that SU(iV c ) gauge theories possess a smooth 
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N c — *■ oo limit if one varies the coupling so as to keep g 2 N c fixed. Moreover the leading cor- 
rection to this limit should be 0(1/N 2 ). It is obviously of interest to check these expectations 
in a fully non-perturbative calculation. In addition, the phenomenology of large- N c theories 
assumes that the theory remains confining in that limit. It is important to check that this 
is really the case. Finally we wish to see how large are the corrections at, say, N c = 3. Is it 
really the case that SU(3) is close to SU(oo)? 

In addition we aim to calculate the detailed mass spectrum of the SU(oo) theory. We note 
that models and theoretical approaches are usually simpler in that limit. For example, the flux 
tube model of glueballs |j, [T| would naively appear to be identical for N c > 2. However, 
because the model does not incorporate the effects of glueball decay, it should really be tested 
against the N c — > oo spectrum since it is only in that limit that there are no decays. A second 
example is provided by the recent progress in calculating the large N c mass spectrum using 
light-front quantisation techniques [Q. Thus our large- iV c spectrum can serve as a useful 
testing ground for models and attempts at analytic solutions of the theory |12, Of course 
in the process we shall also calculate the spectra of the theories at finite N c and these too can 
be used as a testing ground for models and analytic approaches. We remark that the most 
recent examples of the latter [|TJ], |15|, [16|] are intriguingly successful. For example, in [I5"| the 
SU(2) and SU(3) ++ glueball masses are within 15% of our values and the is even closer. 
In [|16j the string tension is calculated for all N c and is within 2% of the values in this paper. 

We now briefly outline the contents of this paper. In the next section we discuss some 
general properties of SU(iV c ) theories in 2+1 dimensions. The aim is not only to set the 
framework for the subsequent calculations, but also to convince the reader that these theories 
are sufficiently similar to their D = 3 + 1 counterparts that what we learn in this paper about 
the former probably tells us something about the latter. We then move on to discuss the 
technical details of how we carry out our calculations. This includes details of our operators 
and of the variational principle that underlies our extraction of excited states from the matrix 
of correlation functions. We also discuss the problem of identifying the continuum spin of a 
particle on a lattice with only cubic symmetry, and provide a simple criterion for doing so. 
In the following section we present our calculations of the string tension. It is here that we 
aim to demonstrate that the theory remains confining as N c — > oo thus providing the crucial 
ingredient for extracting the phenomenology of those theories. It turns out that the string 
tension is the physical quantity that we calculate the most accurately. So it is in this section 
that we shall test most precisely the expectation that there is a smooth limit reached by 
keeping g 2 N c constant, and that the leading corrections to that limit are 0(1/N 2 ). The next 
section contains our calculations of the glueball spectrum. (Since we have only gluons in our 
theories, all the physical particles are colour singlet composites of gluons i.e. glueballs.) We 
shall briefly comment on the features of this spectrum; in particular upon how it compares 
to what we know about the D = 3 + 1 spectrum. It would of course be nice to be able to 
spot some striking regularities in the mass spectrum of the SU(oo) theory. However one can 
only read significance into the details of the spectrum if one has some framework that relates 
those details to the underlying dynamics. For such an analysis, within the context of a model 
in which glueballs are composed of closed loops of chromoelectric flux, we refer the reader to, 



10, 11]. We shall find that the glueball spectra are in fact very similar to each other for all 
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N c > 2. To provide a contrast it is amusing to compare this to the spectrum that one obtains 
in the D = 2 + 1 U(l) theory - the theory that is furthest from U(oo) in the sequence of U(N C ) 
theories. This we do in an Appendix E. We conclude with a brief discussion of our results. 

We have relegated to the Appendices details of some self-contained aspects of our calcu- 
lations which should be of general interest. Appendix A contains a detailed evaluation of the 
effects of over-relaxation in the Monte Carlo update. This is done by directly comparing the 
statistical errors on the masses that one obtains in (realistically large) calculations that contain 
various ratios of over-relaxation to heatbath. We are not aware of any previous comparison of 
this kind. We find that over-relaxation brings a modest but worthwhile improvement, particu- 
larly for the smaller values of the lattice spacing. Appendix B contains some analyses relevant 
to our choice of operators. In particular we study how sensitive is our operator basis to vari- 
ations in the parameter governing the 'blocking'; whether anything is gained by the inclusion 
of operators incorporating 'baryonic' vertices (for N c > 3); and we give more details about 
our variational calculations of the excited states. These are no doubt known to some experts, 
but we hope they will be of use to others in the field. In Appendix C we compare continuum 
extrapolations of the string tension in the bare and mean-field improved couplings. We show 
that using the latter considerably improves the accuracy of the final results. This is of some 
interest because we are not aware of any such previous comparison. Appendix D contains 
details of the calculations we have performed with asymmetric lattices which have timelike 
lattice spacings that are much smaller than the spacelike spacings. The primary purpose of 
this calculation is to check that our procedure for estimating the masses of the heavier states 
is in fact reliable. However our methods for determining the ratio of the lattice spacings are 
likely to be of more general interest and so we develop this calculation in some detail. This 
enables us to calculate the SU(2) spectrum of the theory close to its 'Hamiltonian' limit, and 
to compare it with the symmetric test of universality. Finally, in Appendix E we 

summarise some properties of the U(l) mass spectrum. This is to provide a contrast to the 
SU(iV c ) spectra that are the subject of this paper. 

During the course of this work, we have published some preliminary summaries of some 
of the topics in this paper, as well as on related topics not covered herein. Our early results 
on the SU(2) spectrum appeared in |T7|. (See |nj for an interesting comparison with the 
spectrum of the gauge Ising model.) The preliminary SU(2) string tension was discussed 
in |19| . That paper also contained a study of the width and vibrational properties of the 
corresponding flux tube, which is not repeated here. Also not covered in this paper are the 
SU(2) and SU(3) deconfining temperatures |2(|. (See f2lj for a more extensive discussion.) 
Some preliminary results on the N c dependence in 2+1 dimensions have appeared in ||. This 
paper also contains some very preliminary calculations for gauge theories in 3+1 dimensions. 



2 Some features of D=2+l gauge theories 

In the first part of this section we discuss some of the fundamental dynamical properties 
of SU(iV c ) gauge theories in 2+1 dimensions. Since we are much more familiar with the 
corresponding properties of the same theories in 3+1 dimensions, it will be illuminating to 
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compare the two theories as we go along. 

The second part of the section focuses on the consequences of the fact that in D=2+l 
parity and angular momentum do not commute. We show that this leads to some parity 
doubling in the spectrum; but that precisely what gets doubled is sensitive to the ultraviolet 
and infrared cut-offs that we impose in our calculations. 



2.1 D=2+l - D=3+l? 

Gauge theories in 2+1 dimensions possess a dimensionful coupling: g 2 has dimensions of 
mass and so provides a scale even for the classical theory. By contrast in D = 3 + 1 g 2 is 
dimensionless; the theory is classically scale- invariant. In addition the Coulomb interaction 
in D = 2 + 1 is oc g 2 logr; so the theory is already confining at the 'classical' level - albeit 
only logarithmically. Nonetheless, these apparently quite striking differences are misleading: 
the D = 2 + 1 theory shares with its D = 3 + 1 counterpart its most important dynamical 
properties, as we now remind the reader. 

• Ultraviolet freedom. Both theories become free at short distances. In 3 dimensions the 
coupling, g 2 , has dimensions of mass so that the effective dimensionless expansion parameter 
on a scale / will be 

g 2 3 (l) = lg 2 ^>0 (1) 
In 4 dimensions the coupling is dimensionless and runs in a way we are all familar with: 

In both cases the interactions vanish as / — > 0, although they do so much faster in the super- 
renormalisable D = 2 + 1 case than in the merely asymptotically free D = 3 + 1 case. 

• Infrared slavery. The counterpart of the couplings becoming weak at short distances is that 
they become strong at large distances - "infrared slavery" . This is immediate if we let I | in 
eqns(|l],|2]). Thus in both 3 and 4 dimensions the interesting physics is nonperturbative. 

• The coupling and the mass scale. In 3 dimensions the coupling has dimensions of mass and 
so explicitly sets the mass scale for the theory: 

rrii = Cig 2 (3) 

where mj is any dynamically generated mass in the theory. (For example, a glueball mass.) 
In 4 dimensions the coupling is dimensionless and so, naively, things appear quite different. 
However in fact here too the coupling sets the overall mass scale. It does so through the 
phenomenon of dimensional transmutation: the classical scale invariance is anomalous, the 
coupling runs and this introduces a mass scale through the rate at which it runs: 

mi = CiA (4) 

where A is as in eqn(§). So in both 3 and 4 dimensions the value of the coupling determines 
the overall mass scale. 
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• Confinement. Both theories confine with a linear potential. This is not something that 
we can prove of course. However lattice simulations provide convincing evidence that this is 
indeed the case. Note that although the D = 2 + 1 Coulomb potential is already confining, 
this is a weak logarithmic confinement, Vc(r) ~ g 2 ln(r), which has nothing to do with the 
nonperturbative linear potential, V(r) ~ err, that one finds at large r. 

In addition to these theoretical similarities, we shall see that the calculated mass spectra 
also show some striking similarities. For example: the lightest glueball is the scalar ++ with 
a similar mass, m ++ ~ 4-y/c, in both cases. 

While the above comparisons provide some support for the argument that what we learn 
about D = 2 + 1 gauge theories might have something to teach us about the more interesting 
D = 3 + 1 theories as well, it is important to emphasise that the theories do differ in important 
respects and are certainly not the same. For example, there are no instantons in D = 2 + 1 
non-Abelian gauge theories. This would surely matter a great deal if we were to include 
quarks. Another difference is the fact that the rotation group in two space dimensions is 
Abelian. This has some important consequences to which we now turn. 

2.2 Spin and parity doubling 

In two space dimensions rotations commute: the group is Abelian. So states of spin J do not 
come in multiplets in the way that they do in 3+1 dimensions where the rotation group is 
non-Abelian. We shall use (x, y) for the spatial coordinates and 9 for the angle of rotation. 
We can then define a parity transformation, P, by (x, y) — > (x, —y). We note that the angular 
momentum operator, xd y — yd x , flips sign under parity. That is to say, if some state |0) has 
angular momentum j then the state P\(f>) will have angular momentum —j. 

This last fact has an important consequence for the spectrum. Suppose \j) is some state 
of angular momentum j and energy Ej . Consider the two linear combinations 

\j,±) = \j)±P\j). (5) 

If they are both non-null they will form a pair of states that have opposite parity, since we 
easily see that P\j, ±) = ±|j, ±). Moreover these two states will be degenerate, since P 
commutes with the Hamiltonian H, and so we have the phenomenon of parity-doubling. Of 
course, so far the argument could be equally applied to the case of D — 3 + 1. The crucial 
question is whether both combinations are indeed non-null. Now as long as j ^ the states 
\j) and P\j) are orthogonal because they have eigenvalues ±j respectively with respect to J. 
In that case it immediately follows that neither of the linear combinations in eqn(|5|) can be 
null. This argument clearly fails for j = just as it fails for any j in 4 dimensions. Thus 
we conclude that the j ^ states come in degenerate pairs of opposite parity. While for the 
j = states there is no reason to expect parity doubling. 

The above argument assumes the continuum rotation group. Our calculations, on the other 
hand, will be performed upon a square spatial lattice whose explicit symmetries are rotations 
under tt/2. Does this make any difference? Indeed it does. States of angular momenta ±j 
are distinguished by the phases exp{±ij9} that they acquire under a rotation of 9. We note 
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that for j = 2 in particular, these phases are identical if we restrict ourselves to rotations of 
9 = rar/2. Thus on a square lattice there is no reason to expect parity-doubling for j = 2, any 
more than for j = 0. Of course as we reduce the lattice spacing, a, we expect to increasingly 
recover continuum rotational invariance on physical length scales. The extent to which we do 
will be reflected in the extent to which we recover j = 2 parity doubling in our mass spectrum. 

The rotational invariance is not only broken by the square lattice: it is also broken by the 
fact that our space-time is a finite hypertorus. If the lattice is symmetric in the two spatial 
directions (as it usually will be), this once again leaves us with rotations of ir/2. As the volume 
becomes large compared to the physical length scale of the theory, we expect to recover full 
rotational invariance, and re-obtain j = 2 parity-doubling. 

We thus expect to find parity doubling for j ^ states in the D = 2 + 1 theory. However 
some of this parity doubling may be lost to the extent that either the lattice spacing or the 
periodic boundary conditions affect physical length scales. Thus the restoration of parity 
doubling, in particular for j = 2 states, will provide us with direct evidence for the separation 
of the physical length scale from both the ultraviolet and the infrared cut-offs. 



3 Methodology 

We work on a cubic lattice with periodic boundary conditions. The lattice spacing is labelled 
a and the length of the lattice in the /^-direction is L M in lattice units. The field variables 
are SU(iV c ) matrices. They reside on the links of the lattice and are represented by Ui or by 
Un(n), using an obvious notation. The ordered product of the matrices around a plaquette of 
the lattice is represented by U p or by U^ u {n). We use the standard plaquette action: 

S = P^V- - ±rReTrU p ] (6) 

p ^ V c 

and this appears as a weighting factor e~ s in the Euclidean Path Integral. In the continuum 
limit this becomes the usual Yang-Mills action with 

Note the factor of a that is there because g 2 has dimensions of mass; the dimensionless bare 
coupling, being a coupling on the scale a, is just ag 2 . 

We shall perform a few calculations on lattices with different spatial and temporal lattice 
spacings: a s and a t respectively. In that case we use an action 

S = - -^ReTrU Pa } + - -^ReTrU Pt } (8) 

where the spatial and temporal plaquette matrices, U Ps and U Pt , are multiplied by different 
couplings whose values are chosen to reproduce the desired ratio of lattice spacings, a s /a t . 
How this choice is made is described in detail in Appendix D. 

The main technicalites involve the Monte Carlo and the calculation of masses. We treat 
these in turn. 
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3.1 Monte Carlo 



The Monte Carlo consists of a mixture of heat bath and over-relaxation sweeps. We discuss 
these in turn. 

3.1.1 heat bath 

For SU(2) we use the standard Kennedy-Pendleton |22| heat bath algorithm. This is extended 
to higher groups using the Cabibbo-Marinari [p3] algorithm where effectively one updates some 



of the SU(2) subgroups of the SU(N) matrices. 

An important practical question here is how many of these subgroups to update. Clearly 
the more subgroups one updates the faster we will explore phase space. However one does not 
want to carry this past the point of diminishing returns. To determine an appropriate number 
of these subgroups we have chosen a criterion which involves monitoring how efficiently the 



action of the SU(iV c ) fields is reduced by cooling the fields |24], p5| , when the cooling is applied 
through different numbers of SU(2) subgroups. We recall that to cool an SU(2) lattice link we 
simply replace the matrix that is on that link by the matrix which minimises the action. This 



matrix is easy to determine [24, 25] . A link appears in 4 plaquettes and hence its contribution 



to the action can be written as 

6S t = -|rr{£/,E} (9) 

where the matrix E is the sum of the 'staples' enclosing the link I. Each staple is an SU(2) 
matrix and hence E is proportional to an SU(2) matrix. Then it is easy to see that the matrix 
that minimises 6 Si is given by 

Et 

Ul = W\ (10) 
We note that this is just the choice of matrix that the heat bath algorithm makes if we set 
(3 = oo. Once we have applied this procedure to every link of the lattice we have performed 
a cooling sweep. And we can systematically reduce the action by performing a sequence of 
such cooling sweeps. We can extend this to SU(iV c ) fields by using the Cabibbo-Marinari 
algorithm and cooling within the chosen SU(2) subgroups. In this case the algorithm no 
longer exactly minimises the action. Instead the rate at which it reduces the action is a 
measure of how rapidly it moves through phase space. So our procedure is to generate some 
(plausibly) thermalised SU(iV c ) fields, and then to cool these fields using various numbers of 
SU(2) subgroups. An example of this, for the case of SU(5), is shown in Table [I]. We see that 
if we use very few subgroups the decrease in the action is very slow. (Compared, for example, 
to what happens in the SU(2) theory.) As we increase the number of subgroups the action 
decreases more rapidly, indicating that the algorithm explores phase space more efficiently. If 
we were to increase the number further then clearly at some point we would find that it led to 
little further change in the rate of decrease in the action. At this point we would certainly be 
into diminishing returns. We thus try to choose the smallest number of subgroups that will 
reduce the action reasonably fast. We then use these same subgroups in the Monte Carlo. In 
practice we have used 3, 4 and 8 subgroups in the case of SU(3), SU(4) and SU(5) respectively. 
(There is obviously some ambiguity in the precise choice.) 
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3.1.2 over-relaxation 



In addition to heat-bath sweeps one can also use over-relaxation sweeps [26 , 27| , |28]. In SU(2) 



this corresponds to replacing our old link matrix, U a id, by a new link matrix, U new , defined by 

Unew — Tj^^oWT^T' (H) 

where the notation is as in eqn(^). It is easy to see that this change does not alter the action. 
Moreover it can be extended to SU(iV c ) using the Cabibbo-Marinari algorithm. 

An over-relaxation step involves a large change in the link matrix and so it is plausible 



that it will increase the rate at which we traverse our phase space p6| , |27| , p8| . Indeed, in 4 
dimensions, there is evidence that this is so for large Wilson loops |29| . However what we 
are interested in is the calculation of the low-lieing mass spectrum and so what we want to 
know is how over-relaxation affects such a calculation. In Appendix A we present a rather 
detailed study of this both in SU(2) and in SU(3). (This is, we believe, the only study of this 
kind for gauge theories in 3 or 4 dimensions.) We find that for physical quantities, such as 
masses, a suitable mix of over-relaxation and heat bath sweeps decorrelates field configurations 
significantly, although not dramatically, faster than pure heat bath. There is an additional 
gain that arises from the fact that an over-relaxation step is faster then a heat bath step, 
which in any case involves the calculation of all the staples. (This gain is greater in 3 than in 
4 dimensions since there are fewer staples to calculate in the former case.) 

Thus in the calculations of this paper we shall typically choose to make 4 or 5 over- 
relaxation sweeps for each heat bath sweep. 

3.2 Calculating masses 

Our mass calculations are entirely conventional. The starting point is the observation that 

(0^)0(0)) = J2 I (vac\<f>\n) | 2 exp{-E n t} 

n 

^ |(mc|0|O>| 2 exp{-£ o t} (12) 

where |0) is the lightest state that couples to the operator and E is its energy. (We use 
operators that are localised within a single time-slice.) So if we want the mass of the lightest 
colour singlet state with quantum numbers J,P,C we simply construct a p = operator with 
those quantum numbers, calculate its correlation function and then obtain the mass(=£'o) 



using eqn(|12|). If the quantum numbers are trivial, the lightest state might be the vacuum, in 
which case we use vacuum subtracted operators. Of course it might be that for some quantum 
numbers the lightest state is a multi-glueball state. We shall come back to this possibility 
later on, but shall, for convenience, ignore it for now. 

On the lattice t = na and so what we obtain, not surprisingly, is aE n , the energy in 
lattice units. Note that if we are on a lattice with a finite periodic temporal extent, then 
the expression in eqn(|T2"D needs to have an additional term for the propagation around the 
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'back' of the torus. Such a term will always be included in the numerical calculations of this 
paper, although we shall, for simplicity, persist in writing all our expressions as though the 
temporal extent were infinite. We also note that the temporal extent of our lattice, T = aL t , 
will always be chosen large enough for the partition function, Z, to be accurately given by its 
vacuum contribution: Z ~ exp{— E vac T}. Thus the energies we calculate will always be with 
respect to the energy of the vacuum. 

In principle we can obtain from eqn fll2|) any number of excited states as well. In practice, 
however, fitting sums of exponentials to a function is a badly conditioned problem. So one 
needs to develop a more sophisticated strategy, as described later on in this section. 

Again in principle, one can use in eqn(^) any operator with the desired quantum num- 
bers. However, a numerical calculation has finite statistical errors and because the function 
(0t (£) 0(0)) is decreasing roughly exponentially in t it will, at large enough t, disappear into 
the statistical noise. Thus in prectice we need to be able to extract E from eqn(|12D at small 
values of t. This requires the coefficient |(t>ac|0|O)| 2 to be large. That is to say, we need to 
use operators that are close to the wave-functional of the state in question. 

If we want to use good operators, we obviously need some simple way to decide which 
operator is in fact better. We shall use a variational criterion. However before coming to 
that we say something more about the operators we actually use. This splits naturally into a 
discussion of glueball operators and those from which we extract the string tension. 

3.2.1 operators for glueballs 

We are interested in colour singlet operators because, as we shall see, our theories are confining. 
Now, the trace of an ordered product of link matrices around any closed path of the lattice 
is a colour singlet. So we can build our operators out of such loops. Moreover, under charge 
conjugation the trace will go to its complex conjugate: so the real part is C = + and the 
imaginary part is C = — . For N c ^2we can also construct colour singlet operators containing 
'baryonic' vertices. We shall not use such operators in the calculations of this paper, but 
include a discussion of their properties in Appendix B. 

As a simple example, consider the set of spatial plaquettes U xy (x, t) and form the operator 



It is a colour singlet. Moreover it is translation invariant and so has p = 0. (To obtain a 
non-zero momentum we would include a factor of exp{ipx}.) It is C = + because we take 
the real part of the trace. It is obviously invariant under parity and so is P = +. Finally 
it is obviously invariant under the mr/2 rotational symmetry of our lattice: so it has J = 0. 
This operator will therefore project onto states that have J PC = ++ and p — 0. So from its 
correlation function we can, using eqn(|T2|), extract the lightest ++ glueball mass. 

Suppose we now consider the ordered product of link matrices around an arbitrary closed 
curve C that starts and ends at the point (x,t). Call it Uc(x,t). Then the linear combination 




(13) 



X 



0W = EE^"^eTr{^ (e „ )c 




(14) 
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will have J = j , C = + and P — ±. Here the angles being summed over are 9 n = mr/2. R(6) 
is a rotation operator, so that R(6)C is the contour obtained when we rotate C by an angle 
8. Similarly PC is the parity transform of C. In the second term U is conjugated because the 
order around the curve is reversed under parity. If we replace ReTr by ImTr we get C — — . 
Symmetries of the curve C will be reflected in the operator in eqn(|l4]) being null for some 
values of J PC . 

There is obviously an ambiguity in our assignment of J. We use the continuum notation 
because we are interested in the continuum spectrum and we expect that the lattice will 
recover continuum rotational invariance as a — > 0, at least on physical length scales £/a — ► oo. 
If in constructing our operators we limit ourselves to rotations of mc /2 then any continuum 
spin J that gives the same value of exp{ijn^} will couple to this operator. That is to say, 
the state that we label by J=0 actually contains states with J=0,4,8,12,.. and similarly for 
the states we label by J=l and J=2 (which is all we have with rotations of vr/2). A similar 
ambiguity occurs in 4 dimensions. It is usually assumed that in this tower of states it is the 
state with the smallest value of J that has the smallest mass. Thus, in our case, we shall 
claim to calculate the lightest J=0,1 and 2 states. We shall return to this point later. 

It might be useful to indicate the operators Uc that we actually use. Clearly we need only 
specify the curves C. The first set consists of square and rectangular curves. In particular, 
the lxl, 2x23x3 squares and the 1x2, 1x3, 2x3 rectangles. These curves are 
obviously symmetric under parity reflection. Taking into account that parity also conjugates 
the matrix, it is easy to see from eqn(|l4|) that we can only get J ++ and J states. Moreover 
the square loops can give us only J = while linear combinations of the rectangular loops 
can give both J = and J = 2. However all the loops are symmetric under rotations of it 
and so cannot give J = 1. To obtain J = 1 and P = —C states we need other operators; 
in particular we need curves that are not symmetric under P. To describe such curves it is 
convenient to use an obvious shorthand notation in which the plaquette in the x, y plane 
would be written as xyx^y^. In this notation the curves we use are a path ordered product 
of 2 plaquettes i.e. xyx^y^x^y^xy, and the 'twisted' version of this xyx^y^y^x^yx; the path 
ordered product of the 1x2 loop and a plaquette i.e. xyyx^y^y^x^y^xy, and the twisted 
version of this xyyx^y^y^y^x^yx; and finally the path ordered product of two 1x2 loops 
i.e.xyyx^ V V xxy^ x* x< y and the twisted version of this, xyyx^y^y^y^xxyx^x^ . From suitable 
linear combinations of rotations, parity inversions and real or imaginary parts of these loops 
we can construct operators with J — 0, 1, 2, P — ± and C — ±. 

At this point we have described in some detail the symmetry properties that the operators 
need to have. However all the operators we have described so far are ultraviolet: they are 
based on loops of size 0(a). Such operators will have an approximately equal projection onto 
all states of the specified quantum numbers. The number of excited states increases rapidly as 
a — ► 0. Thus the normalised projection onto the ground state decreases rapidly. This means 
that as a — > we have to go to much larger t in eqn(|T2"D to see the ground state dominating 
the correlation function. But we cannot do so because of the statistical noise in our Monte 
Carlo calculation. This means that we rapidly lose the ability to calculate ground states as 
we approach the continuum limit. 



An efficient remedy for this has been known for a long time [30, 31, 32]. What one needs 
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are operators that extend over physical length scales and are smooth on such scales. Only such 
an operator has a chance of looking like the ground state wave-functional if, as one expects, 
the latter is smooth on physical length scales. Guided by an intuition developed in the context 
of qq wave-functions, one would expect the first excited state to have a node. This could be 
approximated by a linear combination of large smooth operators. Higher excited states would 
be characterised by more nodes. Hence by more complicated linear combinations. While this 
argument is by itself no more than plausible, it turns out that this strategy works remarkably 
well. 

We use the iterative 'blocking' or 'fuzzing' algorithm that has been used extensively in 
D = 3 + 1 spectrum calculations [[H], |32|. We shall not repeat the details here; for a recent 
detailed account (in the context of SU(2) gauge fields coupled to fundamental scalars in 
D — 2 + 1) see [03]. Briefly, at the first 'blocking' level one has the usual link matrices: 
U*(x,t) = U^Xjt). At the second level we construct a 'blocked' link matrix, e.g. U%(x,t), by 
summing the paths xx, yxxy^ and y^xxy and projecting back to the 'nearest' SU(A r c ) matrix. 
All the paths start from the point (x, t) and end at the point 2a away in the x direction. But 
these blocked links are not just longer; they are fatter (in the spatial directions) as well. One 
iterates this procedure: the blocked link matrices are formed in exactly this way from 
the U^" 1 . (All this for spatial ji only.) Thus these operators join sites that are 2 N ~ l a apart, 
and are correspondingly fat as well. We can form path ordered products of these blocked 
links: for example around a super-plaquette, C = xyx*y\ where now each step is of length 
2 N ~ 1 a. The trace of this will be a colour singlet. (After taking expectation values; there may 
be small non-gauge invariant pieces that depend on how we projected from the sum of paths 
back into SU(A' r c ). See Appendix B for a brief discussion.) Clearly the blocking algorithm is 
far from unique. In Appendix B we compare a particular subset of such algorithms in order 
to motivate the particular version we have used. 

Thus we can form large smooth operators on any size-scale we like. When we reduce a by 
a factor of 2, we need only iterate the blocking procedure one extra time. We form operators, 
using a sufficient range of blockings (as determined by preliminary test calculations), on all 
the paths described earlier in this subsection. Thus we often have O(50) different operators 
for any given quantum numbers. Of course we do not need to consider all of these; many are 
dominated by uninteresting ultraviolet excitations. How to choose the 'best' is the question 
we shall return to, after a brief detour describing the slightly different problem of extracting 
the string tension. 

3.2.2 operators for the string tension 

We can calculate the string tension by calculating the energy of the lightest state composed 
of a static q and q a distance R apart. (Any fundamental charges will do; we use q for 
quarks because they are so familiar.) If we have linear confinement then this energy, E m i n (R), 
provides our definition of the string tension, a, as well as providing us with a definition of the 
static quark "potential", V q q(R), 

E min (R) = V qq {R) R ~°° aR. (15) 
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For large R one thinks of this state as being composed of the dressed static quarks with a 
confining flux tube of length ~ R joining them. 

We note that the usual potential that enters phenomenological discussions of the string 
tension |34j] is essentially based on the Schrodinger equation and the relationship with our 



definition is not a simple one; this is apparent if one considers, for example, the case of 



QCD. Vacuum quark fluctuations break the string, so the potential as defined in eqn(15) will 
flatten off for larger R. The phenomenological potential, on the other hand, continues to rise, 
although it may acquire a modest imaginary part to incorporate the decay of the confining 
flux tube. Effectively it incorporates information about the time-scales associated with the 
different dynamical processes that contribute. The two definitions differ most dramatically in 
the large- N c , narrow-width limit of QCD. (Less so in the pure gauge theories of interest here; 
unless, for example, one considers the potential between adjoint sources.) 
To project onto this qq state we define the gauge-invariant operator 

<P(t)=q(0)Y[U iq (R) (16) 

where we can suppose that the quarks are separated along the x-direction and the the product 
of link matrices is along the shortest path joining them. The correlation function of this 
operator, taken from t — to t — T, will, for large enough T, be cx exp{— E min (R)T} . This 
correlation function involves two quark propagators; one from (x = R,t = 0) to (x = R,t = T) 
and the other from (x — 0,t — T) to (x — 0, t — 0). In the m q — > oo limit (which is how 
one implements static quarks dynamically) the quark hops along the shortest available route: 
that is to say its propagator is equal to the product of links along the straight line joining 
its end-points. Thus the correlation function is equal (up to some irrelevant factor) to the 



expectation value of the Wilson loop, (W(R, T)). If we have linear confinement, as in eqn(|l5|), 
then (W(R, T)) oc exp{— E min (R)T} oc exp{aRT}, the usual confining area decay of Wilson 
loops. We can improve this calculation, just as we have improved the glueball calculation, 
using smeared link matrices in eqn([L6|). The timelike link matrices will, of course, not be 
smeared; they arise from the quark propagator calculation. 

We shall use a modified version of the above that employs Polyakov loops rather than 
Wilson loops. Construct a product of link matrices that closes on itself through a spatial 
boundary; for example 

L 

<j) P (x,t) =TrY[U y (x,y + ny,t) (17) 

n=l 

on a L x L spatial lattice. This non-contractible loop is what one gets if one stretches our 
operator in eqn(|I^) till the q and q meet and annihilate. It couples to the corresponding 
state: a flux tube of length L that encircles the torus. Such an operator has zero overlap 
onto any contractible loop. One can readily prove this using the symmetry of the action and 
measure under the transformation U y (x,y ,t) — > z N U y (x,yo,t),Vt where y is an arbitrarily 
chosen value of y and z^- is a non-trivial element of the centre. A contractible loop is obviously 
invariant under this symmetry while the Polyakov loop is not. This argument breaks down 
if the symmetry is spontaneously broken; which occurs, for example, in the high temperature 
deconfining phase. 
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If we sum over x to make (f) translation invariant (p = 0) and form the correlation function, 
we obtain at large t the mass, mp(L), of the lightest state containing a periodic flux loop of 
length aL 

(4(*)Mo)> 



'oc e ~ mp{L)t 



e 



Here we have explicitly included the first correction term which is the translation to Polyakov 
loops |36| of the usual Luscher correction [[37]] for Wilson loops. This correction is 'universal', 



but obviously one needs to test whether the physical flux tube does indeed fall into this 
particular universality class. 

As we have seen above, using Wilson loops produces a heavy-quark potential. This contains 
a Coulomb term which is long range oc g 2 logr in D — 2 + 1. This term will of course be 
screened, but having to disentangle it from the linear piece, at the intermediate values of 
r where the calculations are accurate, can decrease the accuracy of the estimate of a. In 
D = 3 + 1 the Coulomb term is oc 1/r and its presence makes it difficult to identify the 7r/12r 
universal string correction. By contrast, in using as we do correlators of p = sums of spatial 
Polyakov loops, we have completely dispensed with any charges and have transformed the 
problem into a standard mass calculation. Because there are no charges, there is no longer a 
Coulomb contribution. This benefit has of course been achieved at a price: we no longer have 
a calculation of the heavy quark potential, but only of the string tension. 

Just as for glueballs, the simplest operator is too ultraviolet to be useful as a — > 0. To 
remedy this we replace the product of elementary links in eqn(|l7]) with a product of blocked 
link matrices, as defined earlier in this Section. As we shall see, there is always a blocking 
level for which this smeared Polyakov loop is very close to the wave-functional of the ground 
state of a flux tube that winds around the torus. 

Two technical asides. When using link matrices at a blocking level Nb, the sites are 
spaced a distance 2 NB ~ 1 a apart. A given product of blocked links, that starts at say y — 1, is 
not quite invariant under translations in the ^-direction because the blocked links themselves 
are not completely invariant. One can remedy this by summing products that start at y = 
1,2, ..,2 Nb ~ 1 — 1 respectively; and this does in fact improve the operator overlap slightly. A 
second point is that L need not be divisible by the length of the blocked link. In that case 
we include links of a lower blocking level, averaged with staples that include transverse links 
blocked to the level of interest. (In practice this extra smearing with staples is of marginal 
utility in getting a good overlap.) 

3.2.3 variational criterion and excited states 

Our lattice action possesses the positivity properties that allow our lattice correlation functions 
to be decomposed as in eqn flllf ). Let us define an effective mass by: 

(19) 
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Then it is easy to see from the fact that all the coefficients in eqn(|T^) are positive that 



am e ff(t) > am e ff(t + a) Vt (20) 



This is a very useful property; it tells us that m e ff(t) provides an upper bound for the mass, 
ma, of the lightest state with the quantum numbers of the operator 0; whatever the value of t 
and whatever the actual operator used. Since the statistical errors on am e ff(t) increase with 
t, we can assume that any apparent increase of the effective mass with t is in fact a statistical 
fluctuation. 

Now we know from eqn(|l2D that 

am e ff(t) am G (21) 

When is t large enough for this limit to have been effectively reached? Since we know that 
m e ff(t) decreases with increasing t then we can estimate m G by the value of the effective mass 

am G ~ am ef f(t ) (22) 

where t is the lowest value of t for which 

m e ff(to) < m e f f (t > t ) (23) 

within errors. 

Obviously this criterion becomes convincing only if the errors are small enough for the 
relation in eqn([23|) to represent a significant constraint. In practice that will only be the case 
if to is small, which will only happen if we have a 'good' operator; i.e. one which mainly 
projects onto the lightest state. Thus it would be useful to have a simple practical criterion to 
decide, early on in a calculation, which operator is the best. Such a criterion is immediately 
suggested by considering the normalised correlation function: 

rh s = W)0(O)) _ (0 f e- m 0) m , 
U[t) - <0t (0)0(0)) <0t0> W 

Clearly if we were using a complete basis of operators, then the best operator would be the 
one that maximised C(t): it would be the wave-functional of the lightest state and we would 
have C(t) = exp{— met}. If the basis is not complete, this suggests a variational criterion: 
the 'best' operator, 0, is the one which maximises C(t). In practice we shall use t = a. The 
reason is that one obtains an accurate value of C(a) in even a small calculation, and so can 
determine early on which are the operators that one needs to calculate with. The value of 
C(a) provides us with an estimate of exp{-amc} and hence amc, which we know to be an 
upper bound on the true mass. In practice we improve upon this estimate by calculating 
the correlation function of this best operator and getting our mass estimate using the first 
effective mass that satisfies eqn(p3|). 

Our general strategy for obtaining estimates of the ground state and excited state masses 
is an extension of the procedure we have just described. We start with some set of, say, 
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N lattice operators, fa : i — 1,..,N, which we normalise so that (fa* fa) = 1. (These are 
chosen from the operators discusses earlier in this section.) We then find the normalised 
linear combination of the fa that maximises C(a) = (fa(a)faO)). Call this operator $ x . This 
is our best estimate for the ground state wave-functional within the space {(pi}; and the 
associated value of C(a) provides us with a lower bound estimate for exp{— ami} where mi 
is the ground state mass. We can find higher excited states just as simply. First we construct 
a basis of operators, fa : i = 1, .., N — 1, that spans the (N — l)-dimensional subspace of the 
space {<pi} which is orthogonal to $i. We now find the linear combination of these ^ that 
maximises C(a) = (<ft ,J< (a)fa(0)). Call this operator $2- This is our best estimate for the wave- 
functional of the first excited state. The associated value of C(a) provides us with an estimate 
for exp{— am 2 } where m 2 is the mass of the excited state. We can continue this procedure 
obtaining operators <3> 3 , $ 4 , ... from which we can obtain the energies of higher excited states. 

Because our basis is finite the above mass estimates need not be very good. To improve 
upon them we calculate correlations of our approximate wave-functionals, ($J(t)$j(0)}, and 
from these obtain effective masses for as large a range of t as our statistical errors (which grow 
with t) will allow. For each correlation function we look for a 'plateau' in the effective masses 
and use the first mass along that plateau. For the lightest state we are, in principle, looking 
for a plateau that extends to t = oo. For the excited states we expect, with our incomplete 
basis, to have some admixture of lighter eigenstates, and so the initial plateau should be finite 
and will eventually drop to the masses of the lighter states. That is to say, for excited states 
the mass estimate can be lower than the mass of the state whose mass is being estimated. 
This undoubtedly means that there is a larger systematic error on our estimate of the mass of 
an excited state than on that of a ground state. We do not know how to estimate this error 
(for either type of state) but the reader should be aware of its existence. 

We have not yet said how we calculate the $j. We use the following standard procedure 
|. Define the N x N correlation matrix C(t) by 

Cy(t) = (^(^(O)) 

Let the eigenvectors of the matrix C~ 1 (0)C(a) be v 1 ; i = 1, . . . , N. Then 

N N 

*i = Q^ v l k fa = 2 a ik fa (26) 

k=l k=l 

where the constant q is chosen so that is normalised to unity. There are of course many 
variations possible on the above procedure. 

We return now to the choice of our original basis of N operators, fa; I — 1, N. What we 
do is to carry out a short preliminary calculation with typically 5 blocking levels of perhaps 
6 to 12 different operators. We calculate only the diagonal correlation functions. Comparing 
the values at t = a we identify the best operator and a few which are almost as good. We 
also take a number which are significantly worse, since, after all, we want our basis to contain 
a reasonable overlap onto some excited states. The sort of basis that we were easily able to 
accommodate (in terms of memory) had ~ 15 operators. In those cases where we had more 
we split the basis into two and worked with both bases separately. Ideally of course one wants 

15 
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(25) 



to work with a single basis. The smallest basis was for the string tension; but here we were 
only interested in the ground state because by using operators that are translation invariant 
along the Polyakov loop, we automatically exclude any significant overlap onto the interesting 
string excitations of the basic flux loop. 



3.2.4 lattice and continuum J 

Suppose we have an operator <fi obtained by multiplying the (blocked) link matrices around 
some closed curve C. The rotation of this curve by an angle 9 gives the operator <pg. We can 
then form an operator of spin J 

<j){J) = f d6e lje <p B (27) 



in the usual way. This assumes we are in the continuum of course. On our square lattice we 
only use rotations of it/ 2: 

<MJ)=£e lJ ^0 nf . (28) 

n 

As we remarked earlier, 0l( J) is not just spin J but will obviously contain all spins J±4iV, ViV, 
since all these spins provide identical phases at 9 = mr/2. It is nonetheless customary to label 
the lowest energy state by the lowest possible spin, in the expectation that higher spin states 
will naturally be more massive. This is quite unsatisfactory: for example it is really not at all 
obvious that a J=3 glueball must be heavier than a J=l glueball (these are ambiguous since 
J=-l and J=l are degenerate parity transforms and 3-(-l)=4). Which one is heavier can be 
an important issue in any given dynamical model (as, for example, in |l 1]). 

In fact the situation is significantly better than this [ID||, in the case where one uses smeared 



operators with large overlaps onto the ground state. We shall now show this. 

We note that the smeared operators that we construct and, which we then insert into 
eqn(|28|), spread substantially in all spatial directions. We are here only interested in the fact 
that this also involves an angular spread. We might imagine modelling this qualitative feature 
using some function like ~ exp{— 9 2 /a 2 }, with the value of a determining the angular spread 
of the operator. This would be the amplitude to find \9) in (p\vac). The amplitude would 
change to ~ exp{ — (6 — 9 ) 2 /a 2 } if we rotate (ft through an angle 9 . Of course, we cannot 
use precisely this form because it does not reflect the periodic nature of the angular variable. 
However we can modify it slightly so that it does, 

= E e S3_ , ( 29 ) 

n 

and in that case it possesses the Fourier expansion 

-i m=+oo 2 2 

£ e -^e ime . (30) 



m=— oo 



Suppose we now insert this in eqn(^8|) with, for example, J = 0. We obtain 

-i m=+oo 2 2 

MJ = 0) = 77-^ E e-^e ime {l + i m + {-l) m + {-i 
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2 



N=+oo 



E e 




(31) 



N=-oo 



We see, as expected, that we not only have J = but that states with J = AN, ViV also 
contribute. However what is interesting is their overlap, which is oc exp{— 4a 2 N 2 }. We see 
from this that if a is not small, then these higher spin contributions are severely suppressed. 
This should be no surprise: in the extreme limit where our operator is smeared uniformly over 
all angles it is obvious that only J = can contribute. 

We see from this argument that smeared operators will generically have the largest overlap 
onto the lowest |J|. The argument relies on the operator being smooth over some finite angular 
region. This is true of our elementary smeared operators, but is not necessarily true of linear 
combinations of these. Since our variational calculation produces such linear combinations, 
we need to continue the argument a little further. 

What eqn([Jl]) tells us is that states of larger than minimal J will have a suppressed coupling 
to an elementary smeared operator. Thus while it is certainly possible that the lightest state 
with "J = 0" actually possesses J = 4 and that it has a large overlap onto the variationally 
selected linear combination of elementary smeared operators, its overlap onto any individual 
smeared operator should be visibly suppressed. In practice we have found this not to be the 
case in any of the channels: typically we can find a smeared operator for which the overlap is 
> 80%. Thus we can confidently state that the lightest states with J = 0, 1, 2 do indeed have 
those spins, for all values of P, C . We have not attempted to perform a similar check for the 
excited states in these channels. 

Clearly one should use the approximate rotational invariance on scales £ ^> a to construct 
operators that, to a good approximation, have any value of J that one desires. Such a 
calculation is in progress [jSSfl . 

4 Confinement and the string tension 

In the previous section we described how to calculate the mass, mp(L), of a flux tube that 
winds around our L x L spatial torus. Whether such a flux tube actually exists, that is to 
say whether we have linear confinement, will be revealed by how mp(L) varies with L. This 
is the first question we address. 

Having shown that we do have linear confinement, we turn to the problem of extracting 
continuum values of the string tension in units of the mass scale provided by g 2 . This turns 
out to be much less ambiguous than the corresponding D = 3 + 1 calculations where the 
scale is provided by, say, A mom . Nonetheless we shall see that using 'improved' couplings does 
indeed enable us to produce more accurate extrapolations. An explicit demonstration of the 
extent of the improvement is provided in Appendix C. 

Having obtained the continuum values of \/cr/g 2 for the SU(2), SU(3), SU(4), SU(5) gauge 
theories, we then test certain expectations concerning the large- N c limit: 

• is SU (oo) confining? 

• is the N c — > oo limit reached by varying g 2 oc 1/N C 7 
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• is the leading correction 0(1/N 2 )! 

We also get our first indication of how small we can make N c and still be close to the N c = oo 
limit. 



4.1 Testing for linear confinement 



When we are using eqn(|l8D to extract rrip(L), it is not t that we know but n t where t = an t : 



so what we actually extract is amp(L). If we have linear confinement with a string tension a 
then we should find 

amp(L) = a 2 aL + ... (32) 

6L 

for large enough L. Here we have also included the "universal" string correction. Its presence 
is also something we would like to test. 

Since the numerical calculations are fastest in SU(2) that is where we have performed 
our most detailed tests. In Fig.|I] and Fig.|2| we show how amp(L) varies with L for (3 = 6.0 
and (3 = 9.0 respectively. The first thing we note is that there is indeed an approximate 
linear dependence of amp(L) on L, with an apparent trend towards exact linearity at large 
L. Ideally we would like to see this rise continue to L = oo. This is not possible to test in 
a numerical calculation, but what we can ask is whether the linear rise extends to physically 
large values of the string length, aL, or not. Now a convenient physical length scale is given 
by £ s = l/\/<y where we can get a 2 a from the asymptotic linear rise. Doing so we find that 
our largest lengths correspond to aL ~ 8£ s and aL ~ 5^ s at (3 = 6 and (3 = 9 respectively. 
These, we claim, are large distances. For example, they would correspond to ~ 4fm and 
2.5/m respectively in the real world where £ s = \ j\Jo — 0.5 fm. 

We can see from Fig.|l| and Fig|| that the dependence of the mass on L is not exactly 
linear; indeed in the latter figure we plot the ratio amp(L)/L precisely in order to expose the 
deviations from linearity. We note that the approach is from below, i.e. the leading correction 
must have a negative sign, just as it does in eqn(^). We have plotted a fit of this form in 
Fig|| and we see that it appears to be compatible with the observed variation. Can we test 
the correction term in eqn(|32"D more precisely? Suppose the lattice sizes Li are ordered so that 
L i+1 > Li. Let us parameterise the corresponding loop masses by am P (Li) = a 2 a Li — c e ff/Li 
and the same for amp(L i+1 ). Then we obtain 

am P (L i+ i) am P (Li) 

Ceff = (33) 

What are we looking for? At small L higher order corrections in 1/L will be important and 
so c e ff will vary as we increase L. If however c e // — > c as L — > oo then this tells us that the 
functional form of the leading correction is indeed c/L. If the value of c is compatible with 
7r/6 ~ 0.52 then we have some evidence that the correction is of the universal form. How 
much evidence depends on the precision of the comparison of course. 

In Table [5] we list the values of c e ff for various ranges of L iy L i+ i. We also show a single 
value obtained at (3 = 12. In comparing the distances at various values of (3, we can use the 
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fact that lim a ^ (3 = 4/ag 2 which tells us that, roughly, a oc 1//3. So L = 32 — 48 at (3 = 12 
corresponds roughly to L = 24 — 36 at j3 = 9 and to L = 16 — 24 at j3 = 6. As expected we see 
a strong variation of c e ff at the small values of L where our calculations are most accurate. 
As the length, L, of the flux loop increases its mass also increases and so the relative error on 
c e ff increases quite rapidly. So while there is good evidence that for larger L c e ff grows to be 
at least as large as the theoretical value of tt/ '6 ~ 0.52, there is only a little direct evidence, 
from the (3 — 9 values, that this is indeed the asymptotic L — > oo value. Taken as a whole, we 
read the results in Table || as providing significant support for the applicability of the Luscher 
universal string correction to the confining flux tube. We remark that in contrast to D=3+l 
studies using Wilson loops, the present analysis has the advantage of there being no confusion 
with a Coulomb term of the same functional form as the Luscher term. 

Since this is our first serious mass calculation in this paper, it might be worth discussing 
the extraction of those masses in a little more detail. By way of example we list in Table [llj 
the effective masses, as defined in eqn(|l9|), for the /5 = 9, 12 and 14.5 calculations. (The 
last corresponds to our smallest lattice spacing.) We show not only the masses obtained 
using p = operators, but also those obtained using operators with the lowest non-zero 
momentum, ap = 2ir/L on an L x L spatial lattice. From the latter we obtain effective 
energies, aE(p), which we have translated into effective masses using the continuum dispersion 
relation m 2 = E 2 —p 2 . As L decreases, ap = 2it/L becomes larger and at some point it should 
become sensitive to the cut-off at which point this relation will break down. Comparing the 
two sets of masses in Table 10, we observe that they are compatible, within small errors, thus 
demonstrating the restoration of continuum Lorentz invariance. The dispersion relation does 
break down on the L = 6 lattice, but at this point p = 27r/3a ~ 2/a which is certainly an 
ultraviolet momentum. 

It might seem remarkable how small we can make L while still retaining all the string-like 
properties of the flux tube, which after all will have a width of the order of ~ 1/a^/a. In fact, 



as we have argued elsewhere ||19H , this is not surprising if the fluctuations of the tube are not 
too rough, and the transverse volume is periodic. 

We return to the masses. Our criterion is that we choose m(to) as our mass estimate 
if, within errors, m(t ) = m(t)Vt > to- 111 most of the cases shown in Table ITU] that is 
straightforward; the choice of m(t = 2a) seems appropriate. In some cases there is a downward 
drift in the value of m(t) at larger t. For example on the L = 48 lattice at /3 — 12. In this 
case however the mass from the p ^ operator shows no such effect: indeed it shows a slight 
upward drift. This suggests that this drift is a statistical fluctuation. Indeed the large drop 
seen in going from m(3a) to m(4a) cannot be accommodated in any realistic decomposition 
of the correlator that respects positivity. In the case of the L = 32 masses at the same f3 
there is a marginal hint that m(2a) is not asymptotic, but it is difficult to see why it should 
not be if t = 2a is asymptotic at (3 = 9 and at (3 — 14.5 (as it appears to be). The fact 
that correlators can drift [^2J and indeed oscillate [3~9~| outside their apparent errors introduces 
some subjective bias into our analysis. But, as we have seen, it can often be resolved either 
by performing simultaneous p^O calculations, or by considering other calculations at nearby 
values of (3 and L. This renders the problem a minor one in practice, although it may well 
induce a systematic bias at the level of the statistical errors. For this reason we do not take 
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seriously the fact that the mean values of c e ff at larger L tend to be above it/ 6 albeit within 
statistical errors. 

We turn now to the SU(3), SU(4) and SU(5) theories. In each case we have performed an 
explicit test of linear confinement at one value of (3 at least. These are listed in Table |ll[ We 
see that in every case there is an approximately linear increase of the mass with the length of 
the loop. That is to say, we have linear confinement. We also extract and list the values of 
c e ff as defined in eqn([3^). We see that we have some evidence for the validity of the string 
correction in eqn(^) for all our values of N c . 

On the basis of this evidence we shall assume that we have linear confinement at all other 
values of (3 so that we need calculate mp(L) for only one value of L and can then use eqn(|32"D 
to extract a 2 a. Of course aL has to be chosen large enough for the leading correction to be 
the dominant one. We shall use lattices that are about as large as L = 24 at (3 = 9.0 in SU(2). 
As we from Table ||, this should certainly be large enough. 

4.2 Extrapolating to the continuum limit 

In Tables [1^, [1^ and [I4| we list the values of the string tension that we shall use. These have 
been obtained from the calculated flux loop masses using eqn (|32"D . The lengths of these loops 
are also shown in the tables. 

The flux loop masses have been obtained from the cross-correlation matrix, as described 
earlier. The exceptions are all in SU(2): the (3 = 6.56 and L — 16, (3 — 6.0 calculations and 
those at such strong coupling, (3 < 3.47, that the eigenvalue calculation becomes error-driven 
and breaks down. In these cases we applied the simplified variational calculation where one 
chooses, from the original basis of operators, the single smeared Polyakov loop that maximises 
am e ff(t = a). By comparing how such a procedure differs from the full one at neighbouring 
values of L and (3 we believe that any bias induced is within the statistical errors. Once we 
have chosen the 'best' operator, we extract a from am e ff(t = 2a). There are a few cases 
where the naive application of our 'effective mass plateau' criterion would lead us to use 
am e ff(t = 3a) (or larger t). However these are typically two standard deviation effects that 
occur infrequently enough that they can be fluctuations. And in practice if we were to use 
them it would make no material difference to the calculations we now describe. 

We now wish to use these values to obtain the continuum string tension. Since the only 
explicit mass scale is provided by g 2 , we expect that y/a should be some multiple of it. We 
can obtain this ratio from our tabulated values of ayj 'a: 

lim pay/a = 2N C ^- (34) 

(3->oc g 2 

using eqn([7|). 

To perform the limit in eqn (|34D we can add a correction term 

Pay/a = c + ^ (35) 
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and fit the unknown constants, Cq = 2N C } ^- and ci, to the values of Pay/a that we obtain from 



'~9 r 

our Tables. In practice higher order corrections will be important at small f3 and so we will 
need to systematically drop off the lowest-/? points until we get a fit with an acceptable \ 2 - 
Although this is a workable approach, we recall, from the D = 3 + 1 case, that the lattice bare 
coupling provides a poor definition of a running coupling. The basic problems are similar in 
2+1 and 3+1 dimensions and so we might expect that the higher order corrections to eqn(p5|) 
will be much larger than if we were to use a physically motivated coupling. A very simple 
such coupling |4(J is the mean field improved coupling 

(3i = (3x{±TrU p ). (36) 

To define /3j we need the values of the average plaquettes, (4j-TrU p ). These are provided in 
Tables [15], [16] and [H]. In Appendix C we compare extrapolations in (3 and Pi in the cases 
of SU(2) and SU(3), where we have calculations over a wide range of (3 values. We are able 
to demonstrate that the mean field improved coupling does indeed provide a much better 
expansion parameter. Thus we shall extrapolate to the continuum limit using 

Piay/a = c + ^- (37) 
Pi 

in all cases. 



The results of these extrapolations are listed in Table 18 together with the confidence levels 



of the fits and the fitted range. Having obtained the continuum string tensions 

(3* 



0.3353(18) SU(2) 

y/o _ I 0.5530(20) SU(3) 

gi - | 0.7581(40) SU(4) 

0.9657(54) SU(5) 



we turn now to an analysis of their N c dependence. 
4.3 Confinement at large N 

In Fig.[7| we plot our calculated values of ^cr/g 2 against iV c . We immediately observe that the 
variation approaches a linear form for larger N c 

^ oc N c . (39) 
T 

and indeed is nearly linear even down to N c = 2. Now, if our SU(N C ) gauge theories are to 
have a smooth N c oo limit, then in that limit they will have some fixed physical mass scale 
that we shall call fi. If this limit is to be confining we must have 

— -> const ; iV c -> oo (40) 
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From eqn(PU|) and eqn(5C) we immediately infer that 



✓ (41) 

We recall that the usual all-order diagrammatic analysis demands that g 2 oc 1/N C for a 
smooth large-iV c limit. Eqn(£y]) embodies precisely that requirement and so provides a fully 
non-perturbative confirmation of those arguments. 

To complete our demonstration that the theory is confining in the N c — > oo limit we need 
to show that lim^-^oo it^g/V 17 * s finite an d non-zero for the lightest glueball masses. That 
this is in fact the case is something that we shall demonstrate in the next Section; for now we 
shall assume it to be so. 

In addition to predicting that g 2 oc 1/N C , the usual diagrammatic analysis also predicts 
that the leading correction should be 0(1/N 2 ). To test this we fit our string tensions with 
the functional form 

wr« + % (42) 

In Fig.|8] we show how the goodness of fit varies with the power a. From this we can infer that 

a = 1.96 ±0.45. (43) 

If we assume, in addition, that the power should be an integer, then only one value is allowed: 
a = 2. Thus we conclude that we have rather strong evidence that the leading correction is 
also in agreement with the usual diagrammatic expectations. 
Fitting our calculated values, we obtain 

^- 04 975(10)-^ ,44) 



g'N c ■ ' N* 

This fit has a good confidence level, ~ 80%. We note that this tells us something interesting: 
we can describe the physics of SU(N C ) gauge theories, all the way down to SU(2), by that 
of the SU(oo) theory supplemented by the leading correction with a modest coefficient. Of 
course, so far we have only shown this for the string tension: in the next Section we shall see 
that this is also the case for the mass spectrum. 

Before moving on from our result for the string tension, it is interesting to ask whether it 
is possible to quantify the potential error associated with keeping only the leading correction 



in eqn(f44|). There is no unique way to do this, of course, but a first step would be to include 
a higher order correction and see what difference it makes. When we do so we obtain the 
following range of fits: 

> - 04 976(22) - °- 121 ' 43 ' - (45) 



We observe that our result for lim^^oo \f&/ g 2 N c is robust under the inclusion of the higher 
order correction. Our calculations constrain the coefficient of this higher-order correction 
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to be small, and the only significant effect from including it is to double the error on the 
extrapolated value of the string tension. 

The discussion so far has concerned the continuum limit, which is of course what we are 
mainly interested in. However the large- N c expectations will also apply to lattice corrections, 
and we can ask if they are fulfilled. What we would expect is that the 0(ag 2 ) correction 
should also be a function of g 2 N c , i.e. 

<W ° b + b iag 2 N c : iV c ^oo. (46) 



ag 2 N c 

In terms of our fit in eqn(|3T|) this implies that the lattice correction, ci, should be given by 

ci = 4N% : N c -> oo. (47) 



We note that our calculated values of ci, as listed in Table [18], are entirely consistent with 
this being the leading large- N c behaviour. Indeed, if we fit these values with a functional form 
Cx = cN a we find a good fit with a = 4.2 ± 0.6. 



5 The Mass Spectrum 

Having seen that all our gauge theories are linearly confining, we infer that the asymptotic 
states are colour singlet and so we can calculate the mass spectrum using the operators 
described earlier on in this paper. 

We shall first indicate the quality of the lattice mass calculations. We then investigate the 
finite volume dependence of these masses so as to establish control over this potential source 
of systematic error. We shall then carry out the extrapolation to the continuum limit. Finally 
we turn to a study of the dependence of the mass spectrum on N c . We finish with a discussion 
of some features of the calculated mass spectrum. 

5.1 Calculating the masses 

We shall focus on the lightest states because the correlations mediated by heavier states 
decrease so rapidly with t that it becomes hard to know whether we have indeed isolated 
the asymptotic exponential decay. Moreover, glueballs that are heavy enough will decay into 
lighter glueballs and this may require more careful analysis. 

Clearly we want to obtain the ground state in each J channel, and in those cases where 
the ground state is light enough we can estimate one or two excited masses as well. So the 
states whose masses we shall calculate are those of the ++ , ++ *, ++ **, , *, **, 
0-+ 0+-, 2++, 2++*, 2"+, 2"+*, 2—, 2—*, 2+", 2+"*, 1++ 1" and 1+" glueballs. 

We shall calculate the masses, as described earlier, for SU(2), SU(3), SU(4) and SU(5) 
gauge groups and, in each case, for a range of (3 values sufficient to allow a continuum extrap- 
olation. In practice this means for most, but not all, of the values of (3 at which we calculated 
the string tension. 
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For any state, the first question must be: how confident are we that we have indeed 
calculated the mass? That is to say, do we have evidence for an effective mass plateau? 

To address this question we analyse, by way of an example, our SU(5) calculation at the 
highest value of f3. We show in Table |1^ the effective masses we obtain there. Since the 
highest (3 corresponds to the smallest a, these calculations are the closest to the continuum 
limit (and the closest to N c = oo) and so are the ones which are the most interesting. We note 
that it is when the value of a is smallest that the correlation functions drop most slowly and 
we can extract effective masses to larger t. At smaller /3, further away from the continuum 
limit, it will be harder to confirm that we are seeing mass plateaux. 

From this Table we infer that a good estimate of the mass is provided by m e //(2a) in 
each case. That is to say, within errors the effective mass is on a plateau for t > 2a. This is 
self-evident in most cases. In some cases, e.g. for the *, one sees a drop in m e ff of over one 
standard deviation when going from t — 2a to t — 3a. However that is to be expected, just 
statistically, given the large number of correlation functions that we consider. Positivity can 
be useful in such cases. If m e ff(a) and m e ff(2a) are sufficiently close, then one can argue that 
it is not possible for m e ff(3a) to be very much lower. At the margins, this allows us to make 
choices about what is, or is not, likely to be a statistical fluctuation. At this level there is some 
subjective element in the analysis, although it should be evident from Table [19| that this will 
not be an important problem in our calculation. To test this we have performed continuum 
extrapolations using m e ff(3a) whenever a blind application of our criterion for identifying 
mass plateaux so dictated. We found it makes no significant difference although the fits are 
often worse. The reason for the latter fact is that we discount any rise in m e ff(t) with t simply 
because we know from positivity that the effective masses must decrease monotonically with 
t. From the statistical point of view, this is a bias in the procedure which undermines the 
statistical analysis. 

It is apparent from Table O that as we go to heavier states, the evidence for effective 
mass plateaux becomes less significant simply because the statistical errors will overwhelm 
the signal at smaller values of t. As we go to smaller /3, and so larger a, this becomes very 
much worse and we will often not have a useful effective mass beyond t = 2a. In these cases 
we simply assume that m e ff(2a) provides a good mass estimate. This is reasonable. If at a 
high value of (3 a particular operator gives us a mass plateau from t = 2a then at a larger 
lattice spacing, e.g. a' = 2a, an operator that is one blocking level down, and hence half the 
size, should surely give us a mass plateau from t = 2a' = 4a. 

While the above argument is plausible, it cannot replace a direct demonstration. This can 
be provided by allowing the spacelike, a s , and timelike, a f , lattice spacings to differ. We then 
choose a t small enough that the correlation functions fall slowly enough over several (temporal) 
lattice spacings for us to obtain several accurate effective masses. This is an old idea that was 
used precisely for this purpose |^T[ in the early days of glueball calculations. More recently 
it has been used very successfully [^2J as part of the action improvement program. Since this 
is a somewhat different type of calculation to the one in this paper, we leave its discussion 
to Appendix D. The reader will find there an explicit demonstration that even for a coarse 
spatial discretisation, using effective masses at distances between a s and 2a s is an accurate 
way to estimate the masses. 
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In summary, we have taken all the masses that we use in the spectrum calculations of this 
paper from m e ff(t = 2a). We have checked that using m e //(t = 3a) in the few cases that are 
ambiguous makes no significant difference. 



5.2 Finite volume effects 

In a theory with a mass gap, m, and on a periodic spatial volume that is L lattice units 
across, the leading finite size corrections to masses are typically 0(e~ camL ) where the constant 



c = 0(1) will depend on the details of the theory being considered |]43 |. Of course this 
correction will only be relevant once aL is significantly larger than the typical hadronic length 
scale, £. In that case amL will be large, since hadron Compton wavelengths are usually <C £, 
and so the correction will be small: usually too small to be observed with the kind of accuracy 
we possess. This means that we cannot expect to obtain a reliable estimate of the coefficient 
of this correction term. This correction is interesting because it is proportional to a triple- 
glueball effective coupling. However, if what we are interested in is controlling finite-volume 
corrections, then the known functional form of this correction has a very useful consequence. 
Essentially it tells us that if we calculate a mass on volumes aLi and aL 2 which are both 
significantly larger than £, and if we find that the change in mass is small when we go from 
L\ to L2, then we can be confident that any mass shift in going from L = L2 to L = 00 will 
be small compared to the observed change in the mass (as long as a(L2 — Li) ~ £). This is 
important: if the leading correction were power-like rather than exponential then this would 
not be true and controlling finite- volume corrections would be appreciably more difficult. 

So our strategy to control finite volume corrections is as follows. We calculate masses on 
a range of lattice volumes. We include volumes that satisfy the conditions of the previous 
paragraph. And once we observe very small changes on our larger volumes, we can be confident 
that the mass calculated on the very largest volume is identical, within errors, to the L = 00 
mass. 

In practice it would be wasteful to perform such an analysis at each value of (3. Instead we 
choose a couple of values of (3 where we perform an extensive analysis, including very large 
volumes in order to make sure there are no unpleasant surprises. This allows us to establish 
what volumes are large enough that any change in mass becomes invisible within the typical 
statistical errors of our calculations. We then use scaling to infer how this translates to other 
values of (3. 

We perform these calculations in SU(2), simply because that consumes much less computer 
time. Having determined how large the volume has to be in, say, units of 1/a^/a we can take 
this criterion over to SU(3) etc. Of course there is some danger in doing this and so we perform 
at least a modest finite-volume check for each of our non-Abelian groups. 

As a first step, we show in Fig.|9| how some of the lightest masses vary with the size, L, of 
the spatial volume, in the case of SU(2) and at a coupling of (3 — 9.0. This, as we see from 
Table 0, corresponds to quite a small lattice spacing. Our spatial length varies from L = 6 
to L = 32 which corresponds to a variation of ~ 1 to ~ 5 in units of 1/a^/a (about 0.5 to 2.5 
fermi if we were in QCD); a range of sizes that satisfies the conditions laid out above. Note 
that since this is SU(2), we have no C = — states; these shall appear in our (less extensive) 
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SU(3) study. So what we show in Fig|5]are the masses of the lightest ++ , 0~ + , 2 ++ , 2~ + , 1 ++ 
and 1 h glueballs. We also show twice the mass of the periodic flux loop, 2amp, for reasons 
that shall soon be apparent. 

There are several observations we can make from Fig.|9| and the calculations on which it is 
based. 

• As we decrease L we do indeed observe the onset of substantial finite size effects. 

• The parity doubling that we see at large volumes is badly broken by these finite-volume 
corrections. This is not unexpected: the toroidal boundary conditions break the effective 
rotational symmetry from the (dynamically restored) continuous one down to rotations of 
7r/2. As discussed earlier, this undermines the argument for 2 ± degeneracy although not 1^ 
degeneracy. This is precisely what we observe in Fig.[5| Thus the observed degeneracy of the 
2 ± states can serve as a criterion for the lack of finite volume effects. 

• We observe that the value of L at which the 2 + begins to show finite- volume corrections is 
roughly where the asymptotic glueball mass equals twice the flux loop mass: 

2am P (L) ~ am 2 +/ L=00 . (48) 

The same is true for the + . Since the latter is lighter, and since mp t as L f, the scalar 
becomes volume-independent at smaller volumes than the tensor. 

This correspondence with 2mp is easy to understand. Suppose we denote by l x the p = 
smeared Polyakov loop in the x-direction which has the best overlap onto the x-periodic flux 
tube. Typically this overlap will be ~ 90 — 100%. Consider now the operator l x l x - This will 
also be p = and colour singlet, but it falls into the sector of contractible loops and so can 
couple to glueball states. If the transverse spatial size, L y , were very large, then this operator 
would mainly couple to a state that consists of two periodic flux loops whose energy would 
be 2mp(L). We shall refer to such states as 'torelons'. On our lattices L x = L y = L and as 
L becomes small these flux loops will necessarily interact; thus the lowest energy will deviate 
somewhat from 2mp. We can form ++ and 2 ++ combinations, l x l x + l y l y and l x l x — l y l y 
repectively. Again, on large volumes these will mainly couple to states with two flux tubes 
and mass 2mp(L). On smaller volumes, aL < £, the interaction between the flux tubes will 
split these states away from each other and from this mass. While we cannot predict the 
precise variation of these masses with L, one would expect them to decrease, at least until 
aL <C £. Thus the fact that the 2 ++ mass begins to decrease with decreasing L just when 
eqn (|48"l) is satisfied, would seem to simply reflect the fact that for smaller L than this the 2 ++ 
state which is composed of a pair of flux loops becomes the lightest state in that sector. And 
the same for the ++ . 

We have explicitly confirmed this scenario. For large L at least one of our usual glueball 
operators has a large, ~ 90 — 100%, overlap onto the lightest ++ or 2 ++ state. By contrast 
the double flux loop operators have poor overlaps. The value of L at which the lightest mass 
begins to decrease as L decreases, marks the point at which things reverse. For smaller L it 
is one of the double flux loop operators that has a very high overlap onto the lightest state 
and the usual operators all become poor. 
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We remark that similar finite volume effects are observed in D=3+f gauge theories flH|; but 



because mc/^/a is slightly higher in D=2+f, the effects occur on somewhat larger volumes, 
and so their interpretation is that much less ambiguous. 

• Naively we would expect the spatial size at which we begin to encounter large finite size 
effects to be related to the size of the glueball. However in the case of the ++ and 2 ++ 
glueballs, we have seen that the onset of finite size effects is simply determined by the (infinite 
volume) mass and the string tension. Thus the fact that we can go to smaller spatial volumes 
for the scalar than for the tensor, before encountering large finite size effects, is not telling us 
that the size of the scalar is less than that of the tensor. The same holds true in the case of 
four dimensions. 

• In addition to the above, there are finite size effects, visible in the 0~,2~, 1 ± states, whose 
onset appears at much smaller values of L, and which does not appear to be linked to mixing 
with torelon states. (There are no simple torelon states with these quantum numbers.) 

• We infer from Fig.|9| that at (3 = 9 a spatial size of L = 24 is large enough for the lightest 
glueballs to be free of finite size effects within our statistical errors. Assuming scaling and 
(3 = 4/ag 2 this implies that at a general value of (3 a safe size is L > 24(3/9. One can do 
better by using the calculated value of a\Ja to set the scale. Doing so, one can then extend 
the criterion to SU(3) etc. 

In practice we are more cautious than this and have performed an extensive finite volume 
analysis at (3 = 6 as well. This is for a larger value of a, and we include larger lattices ; up 



to more than 8 in units of l/a^/a. This masses are shown in Table [2CJ. We also have a more 



limited study, on L = 32 and L = 48 lattices, at (3 — 12 where a is smaller. These masses will 



be displayed later, in Table |4], where we display our U V = oo" SU(2) mass values. All these 



results confirm our criterion for what constitutes a safe volume. 



In Table |21] we show our SU(3) study. The presentation here is slightly different to that 
in Table in that we show estimates of the ++ and 2 ++ torelon masses using the operators 
described earlier on in this Section. The glueball masses have been obtained from the usual 
glueball operators based on contractible loops. (Of course, the torelon and glueball operators 
do mix and at t — > oo we would always find the same effective mass. But if the mixing is 
small one will, in general, find different effective mass plateux at small values of t.) Just as we 
saw in the case of SU(2), it is clear that the onset of large finite size effects for the scalar and 
tensor glueball masses is linked to the mass of the corresponding torelon state. We are also 
now able to see what happens to the C = — states and we show the which is the lightest 
of these. Since our a is not small, we do not have accurate values of m e //(t = 2a) for the 
heavier states. (As usual our quoted masses are extracted at t > 2a.) We can of course look at 
m e ff(t = a) which has the disadvantage of having some excited state component, but which is 
accurately calculated even for the heaviest states. Although we do not show the values here, 
we remark that they show no sign of any finite size effects that violate our above criterion, 
either for the C = + or for the C = — states. The same is true for our modest SU(4) and 



SU(5) finite volume studies, which appear in Table ^2 



This establishes the level of our control of finite volume effects. There is one further 
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important point. In the case of the 2 ++ the torelon appears to exist as a bound state in 
the mass spectrum for larger L. This makes it difficult to extract a consistent picture of the 
excited 2 ++ states. Since there is no such difficulty for the excited 2~ + states, and since these 
states should be degenerate with the non-torelon 2 ++ excitations, we do not try to overcome 
this difficulty. So the reader should not be surprised to find no masses being quoted for the 
excited 2 ++ states later on. As an addendum to this, we remark that this 'difficulty' appears 
to disappear for larger values of N c . We assume that this is a manifestation of the suppression 
of all mixings at large N c . 



5.3 The lattice mass spectra 

In the previous two subsections we discussed our criteria first for minimising the systematic 
error associated with the extraction of masses on a given lattice, and second for controlling 
the finite volume corrections to such masses. We now use those criteria to extract our "infinite 
volume" lattice mass spectra. 



In Tables £3|, |24|, [25|, ^6] and |27| we list some of the masses that we have extracted in the 
SU(2), SU(3), SU(4) and SU(5) calculations respectively. All are in lattice units. In addition, 
all these masses have been extracted from am e ff(t = 2a). There are a few exceptions to this. It 
sometimes occurs, particularly at the smaller values of (3, that am e ff(t = a) < am e ff(t = 2a). 
On the other hand, we know from positivity that am e ff(t) must decrease as t increases. Since 
the error at t = a is smaller than at t = 2a, it is clear that the t = a effective mass is the 
better mass estimate in these cases, and that is the value we list. However the error quoted is 
the larger one which is associated with t = 2a. In principle when we come to extrapolating to 
the continuum limit we should use the t = 2a effective masses in these cases since otherwise 
there is a systematic downward bias in the statistical analysis. (That is: we correct some 
large upward statistical fluctuations in am e ff(t = 2a), but none of the ones that are large and 
downwards.) However so few values are affected, and these are usually at the lowest values of 
(3, that we choose not to complicate the analysis by doing so. We also note that the volumes 
used here are always at least as large as the minimum necessary, as indicated by our earlier 
finite volume studies. 

We begin with a brief technical aside. Our calculations, at each value of (3 and L, typ- 
ically involved between 80000 and 200000 Monte Carlo sweeps with calculations of glueball 
correlators being made every 5 sweeps. Typically we would have 3 or 5 over-relaxed sweeps 
for each heat bath sweep. By comparing the values of a^Jo the reader can see that the SU(4) 
and SU(5) (3 values are more-or-less equivalent. They are also nearly equivalent to some of 
the SU(3) and SU(2) p values. 



Some comments now, starting with the SU(5) masses listed in Table |27]. Here we focus 



on features that might affect the reliability of the calculations; we leave a discussion of the 
physics till later. We first note that the lightest J = 1 and J = 2 states display parity 
doubling within errors. In contrast to the marked lack of doubling in the J = sector. This 
confirms that we have made our ultraviolet cut-off small enough, and our infrared cutoff large 
enough. The same is true of the J = 2 excited states. However it is less clear what is going 
on amongst the excited states in the J = 1 sector; there appears to be a near-degenaracy 
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between the ground and excited states in some cases. And the expected degeneracy between 
the 1++* and the 1 h * appears to be broken. This may indicate the presence of a J = 3 state 
which is nearly degenerate with the J = 1 state; or it may be that some of these states are 
multiglueball scattering states; or there may be finite volume corrections. We are not well 
placed to distinguish amongst these possibilities in our present calculations. For example to 
investigate the last possibility, we need to do a finite volume study for a small value of a where 
these very heavy masses can be accurately calculated. The only calculation of this kind is in 
SU(2) at (3 = 12 (see Table |24]). We do not see any trend for the 1 ++ * and masses to 
converge as we increase L from L = 32 to L = 48. So it does not seem to be a finite volume 
effect. If we compare different /3 values there appears to be no trend for this effect to decrease; 
so it would not seem to be a finite-a effect. This is a puzzle. As far as the 1 ±_ states are 
concerned, we have even less to go on, because we have no C = — in our SU(2) studies. These 
oddities need to be resolved but they afflict the very heaviest of the states we study and so 
we shall not attempt to resolve the issue here. 

Although we have carefully chosen the volumes so that the ground state J = and J = 2 
masses are essentially infinite volume, this is not necessarily the case for the J = 2 excited 
states. Indeed we observe in Table 21 that the mass of the 2 ++ * is volume dependent and 
is only degenerate with the 2 h * on the largest volume, and then only in those cases, (3 = 9 
and 12, where this is exceptionally large. Explicit calculations with the double Polyakov loop 
operators described in the previous subsection indicate that this is an artifact of the presence 
of a corresponding scattering (or bound?) state whose mass increases approximately linearly 
with L. So as L f the mass moves out of the range of masses we probe. We remark that we 
have no clear evidence of a corresponding J = state. (It would, in any case, not interfere 
with the lowest ++ excitations because these are so light.) It is interesting that this problem 
appears to disappear for the larger SU(N) groups. A possibility is that, as expected, the mixing 
between the double flux loops and our ordinary 'local' glueballs is suppressed by powers of 
1/JV. To investigate this properly one needs to include both double flux loops and our normal 
contractible loops within a single basis for our cross-correlation matrix (as has been done in 



45] for precisely this purpose.) We have not been able to do this here because of the very large 



storage costs this would have entailed. Since we cannot resolve the states unambiguously, we 
shall not attempt a continuum extrapolation of the SU(2) 2 ++ *. For larger groups it seems 
that this problem is not there and so we shall attempt to obtain the corresponding continuum 
masses. 



5.4 The continuum mass spectrum 

A lattice spectrum is only interesting insofar as it can lead us to the spectrum of the corre- 
sponding continuum theory. To obtain the continuum spectrum we need to extrapolate our 
lattice masses to a = 0. The first step is to take ratios of masses so that the scale, a, in which 
they are expressed cancels. We choose to take ratios of the glueball masses, amc, to a^/a 
since the string tension is our most accurately calculated quantity. 

The second step is motivated by the observation |^6| that in pure lattice gauge theories 
the leading lattice correction to dimensionless ratios of physical quantities, such as mc/y/cr 
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is 0(a 2 ). So for small enough a we expect the a-dependence to be given, just as in four 
dimensions, by 

+ ca 2 a. (49) 



y/a(a) y/a(a = 0) 

Of course, instead of using the correction term ca 2 a we could use ca 2 rfiQ where arhc is any 
calculated glueball mass. The difference is formally 0(a A ). The reason for choosing a is simply 
that it is so accurately determined. An alternative way of extrapolating such a mass ratio to 
a = is to use the fact that g 2 has dimensions of mass, and that lirn8_ KX) /3 = 2N c /ag 2 . Thus 
for small enough a 

m G {a) _ m G (a = 0) c_ 

Vv(a)~ s/a(a = 0) (3 2 ' 1 j 

What we do in practice is to choose one of the above forms and attempt to fit all the mass 
values of some state with it. If a good fit is not possible we assume that this is because the 
largest value of a used is too large for the 0(a 2 ) correction to be adequate. So we drop the 
mass corresponding to the largest value of a and try again. We keep doing this until we get 
a good fit. 

In Fig.|ll| we show some examples drawn from the SU(3) calculation. Since the mass ratios 
are plotted against a 2 a, the continuum extrapolations, in eqn (|49"D , are simple straight lines. 
The really striking feature of this plot is how little variation with a there is. This will make 
for unambiguous and very accurate continuum extrapolations. 

Our continuum extrapolations for our various theories are displayed in Table EHl They 



have been obtained by fitting the form in eqn(BOf) to the masses listed in Tables E3L 24 



and the string tensions listed in Tables |T^, [T3], [TJ]. The quality of each fit, as given 



by the confidence level, is given in Table |29|. We have also performed extrapolations using 
eqn (|49|) ; these give essentially identical results, with any differences being much smaller than 
our quoted errors. 

The reader will note that the mass of the 2 ++ * is missing from the SU(2) and SU(3) 
columns. This is because we could find no acceptable fits. We believe this is related to mixing 
with torelons, as discussed previously. The mass of the latter depends sensitively on the 
volume and so will not provide a consistent set of masses at different values of (3 (since the 
volumes are not exactly the same). This is no longer a problem with SU(4) and SU(5) and we 
assume that this is because any such mixing becomes suppressed at large N c . A calculation 
including overlaps between torelon and glueball operators would resolve this question, but we 
have not carried this out. Why we cannot get an acceptable fit for the SU(2) 2 h * is less 
clear. The reason might be that our SU(2) calculations in the 2 h channel had a smaller basis 
(4 operators) than in the later calculations with larger groups. This meant a very small basis 
for the excited states. 

In addition to these spectra we have performed in Appendix D some calculations with 
a very asymmetric lattice action, a t — a s /4. This may be thought of as being close to the 
'Hamiltonian' limit, and it is interesting, as a test of universality, to confront this spectrum 
with the SU(2) spectrum that we have obtained in this section. This we do in Table ^ and 
we observe good agreement within errors. For a detailed discussion of our calculations with 
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the asymmetric lattice action we refer the reader to Appendix D. 

We now have all our continuum spectra and can turn to their dependence on N c . 



5.5 The N c dependence of the mass spectrum 

We can already see from Table that the variation of our mass ratios with N c is weak; and 
that it appears to become weaker with increasing N c . 

To illustrate this we plot in Fig.|l2| and Fig.[l3] the quantity mg/ 'g 2 N c , which is obtained 
from the ratios in Table EH] and the string tensions in eqn(|38|). We choose to plot against 1/N 2 



because the usual diagrammatic analysis predicts that at large enough N c we should expect 

00 ~MV ( 51 ) 



g 2 N c 00 N% 

where Roo = limjv^oo fna/g 2 N c . So on our plot this will be a simple straight line. We observe 
that the dependence of our masses on N c is really very weak indeed, all the way down to 
N c = 2. This indicates once again that the mass scale of the SU(iV c ) theory is oc g 2 N c , as 
expected from the diagrammatic analysis. 



In Table |33] we list the N c — > 00 limits and the slopes, Ri, that result from fitting our 
continuum masses with eqn(|5l|). We can, of course, perform a similar analysis using ma/ ' \J<J 
instead. The results of the corresponding extrapolations are presented in Table [34]. 

Is there anything we can add to our previous result, in eqn(|3|), on the power of the leading 
correction? The only mass that is accurate enough to be potentially useful is m ++ /g 2 N c . 



However, as we see from Fig.|T2], this varies almost not at all with N c and so provides us with 
no useful information on the power of this correction. (The stronger variation in m ++/^/a 
simply reflects the variation of ^/a which we have already studied.) 

For purposes of comparison, it would be interesting to provide an example of a mass 
spectrum that is quite different to the one we have calculated here and yet comes from a theory 
with a dynamics that is not so dissimilar as to make the comparison meaningless. A natural 
possibility is to consider the U(l) theory. Since the leading-order large-A^ arguments are in 
fact for U(N C ), this theory belongs naturally to the sequence of theories we have considered. 
And yet it is so far from N c — 00 that we would not expect it to fit into the pattern we 
have observed so far. As far as dynamics goes, it is a lattice gauge theory which is linearly 
confining and free at short distances. We discuss our results for the U(l) mass spectrum 
(listed in Table |35|) and some peculiarities of the theory, in Appendix E. Here we merely note 
that in the U(l) theory the ++ is no longer the lightest state; the is about half its mass. 
We also note that the mass ratio m ++/^/ cr is much lower than in SU(2) (or in SU(oo) for 
that matter). So this spectrum is indeed quite different; and the comparison enhances, by 
contrast, our claim that SU(2) ~ SU(oo). 

The lack of any visible N c dependence in our most accurately calculated mass, m ++ /g 2 N c , 
is quite striking and provides strong evidence that there is a smooth nontrivial large- N c limit, 
with a physical mass scale oc g 2 N c . Coupled with our previous analysis of the string tension, 
this also tells us that the N c — 00 limit possesses linear confinement. (As we see immediately 
from the behaviour of rriG/y/cr.) From the intercepts and slopes listed in Tables and [33 
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and from eqn(|3), we can obtain the mass spectrum for any value of N c . In this very concrete 
sense we can say that SU(iV c ) theories are close to SU(oo) all the way down to SU(2). Thus 
the large- iV c analysis unifies our understanding of all SU(iV c ) theories in a compact and elegant 
fashion. 



5.6 Features of the mass spectrum 

The purpose of this paper is to calculate the mass spectrum. Extracting interesting physics 
from the detailed features of that spectrum is something that belongs elsewhere. However it 
would be churlish of us not to make a few comments. These will be brief and incomplete. 

• There is clear evidence for the expected parity doubling in the cases of the 2 ±+ , 2 ± ~, 2 ±_ *, 
1 ±+ and 1 ±_ . In contrast, for the J = states, where we do not expect parity doubling, the 
splitting between the P = + and P = — states is huge. 

• Our lightest glueball state is the ++ ; just as it is in D=3+l gauge theories. Moreover its 
mass, in units of \Ja is not that different. If we take the SU(3) continuum extrapolation in 
p5| , which uses the D=3+l lattice glueball mass calculations in fl32| , f49f , and if we perform a 



(52) 



corresponding SU(2) continuum extrapolation using the D=3+l lattice calculations in |47|, |50| 
then we find 

m 0+ + f 3.87(12) SU(2) ; D=3+l 
~ \ 3.65(11) SU(3) ; D=3+l 

The fact that the D=3+l mass ratio is smaller than the one in D=2+l follows naturally 
PH |TTf in flux tube models of gluonic states ||. (It does so from the fact that the closed flux 
loop has more transverse dimensions in which to oscillate; this increases the corresponding 
"Casimir energy", and so decreases the mass of the loop, for a given loop length.) We also 



note from eqn(p2|) that the D=3+l N c dependence has the same sign as in D=2+l. 

• Just as in D=3+l fl2~5l, 152], HSfl, the next heaviest state in the C = + sector is the 2 ++ 



(ignoring excitations of the ++ since these have not been calculated in 4 dimensions). The 
scalar-tensor mass splitting is not dissimilar: e.g. 

m 2++ j 1.65(3) SU(3) ; D=2+l 

m 0++ 1 1.41(7) SU(3) ; D=3+l [b6) 



• Unlike the C = + sector, the C = — sector is very different from its D=3+l counterpart. 
For example, we have a light , while there are no light C = — states in 4 dimensions. 
This may arise from the fact that in 3 space dimensions there is an interplay between C and 
J that does not exist in 2 space dimensions. Consider, for example, a circular flux string. It 
will have an arrow on it, for N c > 3. Under C the direction of the arrow flips. In 3 (but not 
2) space dimensions we can rotate the circle by 7r around a diameter and this also flips the 
direction of the arrow. Note that this means that a rotationally symmetric linear combination 
of such circular loops cannot be C = — . One needs a fluctuation away from a circle to allow 
J = and C = — and this raises the energy. Of course we have gradually incorporated some 
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dynamical assumptions as we moved through the last few sentences. One needs to make the 
argument within a specific model framework and that belongs elsewhere || [10], |TT|| . 

• We observe that whatever splits the C = + and C = — states is weakly dependent on N c ; 
and survives the N c — > oo limit. On the other hand, in SU(2), where we have no C = — , the 
spectrum is clearly a smooth continuation of the N c > 3 C = + spectrum (since our simple 
mass fit encompasses 2 < iV c < 5). This provides a constraint on dynamical mechanisms for 
the C = ± splitting. 

• There are some striking approximate degeneracies in the spectrum. The typical pattern is: 
mo++« — m — , 7Bo++** — m — * and similarly for the J = 2 states. Again, if this is not an 
accident, it does suggests some simplicity in the dynamics. 

To go further requires confronting specific models with the spectrum we have calculated 
here. That goes well beyond the scope of this paper. 

6 Conclusions 

In this paper we presented our calculations of the mass spectra and string tensions in three 
dimensional SU(2), SU(3), SU(4) and SU(5) lattice gauge theories. From these we obtained 
the corresponding continuum spectra. The accuracy of these continuum results reflects the 
large range in the lattice spacing, a, over which we performed our lattice calculations. We 
can compare this range to that in the more familiar D=4 SU(3) gauge theory by using the 
calculated values of the string tension, a 2 a. Doing so we observe that the (useful) range of our 
D=3 calculations would correspond to 5.50 < (3 < 6.55 for the case of SU(2), 5.50 < (3 < 6.50 
for SU(3), and 5.70 < (3 < 6.35 for both SU(4) and SU(5). This range, and the statistics of 
our calculations, is the primary reason why our D=3 calculations are so much more accurate 
than what is available in four dimensions. We also gain something from the fact that our best 
operators are slightly better in three than in four dimensions. 

We noted some strong similarities, in the C = + sector, between the 2+1 and 3+1 dimen- 
sional spectra. This should provide an interesting test for models of glueballs. Indeed one 
of the main motivations for our calculations is to provide a detailed spectrum against which 
models and analytic approaches can test themselves. 

At the more technical lattice level, we studied, during the course of our calculations, the 
effectiveness of over-relaxation, the use of asymmetric lattice actions, how good are operators 
with baryonic vertices, the efficiency of our 'blocking' algorithm, and the extent to which the 
mean-field/tadpole improvement of the coupling really represents an improvement. In this last 
case we were aided by the super-renormalisability of the theory; this allowed us to compare 
directly extrapolations using the improved and bare lattice couplings, in a way which is not 
at present possible in four dimensions. 

The primary purpose of our calculations was to study the large- N c limit of SU(iV c ) gauge 
theories in 2+1 dimensions, and to compare the results of our fully non-perturbative calcula- 
tions with the standard expectations obtained from all-order perturbation theory. We found 
that there does appear to be a smooth N c — > oo limit and that it is obtained, as expected, 
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by varying g 2 oc 1/N C . The leading correction is 0(1/ N 2 ), again as expected. We found 
that confinement - the crucial ingredient for the usual phenomenology - does indeed survive 
the large- N c limit. And we obtained the detailed mass spectrum in that limit. Finally, we 
observed that even SU(2) is close to SU(oo), in the sense that the difference between the mass 
spectra can be described by just the leading 0(1/ N^) correction. 

Thus all D = 2 + 1 SU(iV c ) gauge theories can be described by the SU(oo) theory with 
a modest 0(1/ N"^) correction. This provides a very elegant way to unify and understand all 
these potentially quite different theories. 

There is a wealth of large- N c expectations that we have not explored. For example those 
involving decays, G — ► GG, and, more generally, the N c - dependence of matrix elements in- 
volving various products of singlet operators, as well as their factorisation properties. Neither 
have we attempted to expose the existence of Witten's Master field or to determine its 
properties. All these topics should be readily accessible in three dimensional calculations of 
the kind presented in this paper. The reason we have not addressed them in this paper is not 
because we find them less interesting than the questions we have addressed, but because the 
SU(2) and SU(3) calculations were completed before we realised that we might have something 
interesting to contribute concerning the large- N c limit. 

These calculations also need to be extended by the inclusion of matter fields in the funda- 
mental representation. In this case the leading corrections are expected to be larger, 0(1/N C ), 
and so it is an interesting open question whether SU(2) or even SU(3) will remain close to 
SU(oo). Needless to say, all the above questions need to be addressed in four dimensions; and 
what we can say there is that the first indications j8| are quite promising. 
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Appendix A : Testing the benefits of Monte Carlo over- 
relaxation 

Although there is no systematic procedure known for reducing the exponents associated 
with critical slowing down in D = 4, or D = 3, non-Abelian gauge theories a method that 
appears to have some effectiveness, and which is now in common use, is to mix heatbath and 
overrelaxation J2(], sweeps during the update (for reviews see [p8"|). As far as testing the 



efficiency of this method is concerned, what is available are studies of the decorrelation of 
blocked Wilson loops in SU(2) and SU(3) which show that there is a strong reduction in 
fluctuations when most heatbath sweeps are replaced with overrelaxation. This has helped to 
motivate the widespread use of overrelaxation. 

However, useful as these tests are, what one would like to see is how the statistical errors 
on the physical quantities of interest (glueball masses, string tension, ...) are reduced when 
some fraction of the heatbath sweeps are replaced with overrelaxation sweeps. In D=3+l such 
an exercise would be prohibitively expensive for the small lattice spacings where the answer is 
interesting; and so, as far as we know, no study of this kind has been published. In D — 2 + 1, 
such a study becomes possible and this is what we shall present in this Appendix. Because 
the D=2+l and D=3+l theories have so much in common - in particular they both become 
free at short distances - we can hope that what we find has some relevance to four dimensions 
as well. 

Our heatbath and over- relaxation algorithms have been described in Section 3.1 . We note 
that in both SU(2) and SU(3) the overrelaxation algorithm explores phase space at a constant 
value of the total action. We shall characterise the update pattern by the ratio, R Q , of the 
number of overrelaxation sweeps to the number of heatbath sweeps. Since we use a pipelined 
CPU, all our sweeps employ a variation of a chequerboard update. 

Our study of SU(3) is the more extensive of the two and so this is where we shall begin. 
We have performed comparisons at three values of (3: at — 11 and 15 on 12 2 16 lattices, and 
at j3 = 21 on a 24 3 lattice. If we use the calculated string tension to set the scale of the lattice 
spacing, then these three values of (3 correspond to (3 ~ 5.7, 5.9, 6.15 respectively in the D = 4 
theory with which the reader is probably more familiar. The lattices at (3 — 11 and (3 = 21.0 
are effectively of infinite physical volume for the quantities we shall be considering here. The 
lattice volume at (3 = 15 is of an intermediate size, which mainly effects the nature of the 2 ++ 
glueball. At (3 = 11 we performed calculations for R Q = (pure heatbath) and for R = 5. 
Each calculation involved 80000 sweeps with the data split into 40 bins of 2000 sweeps each 
for the error analysis. At (3 = 15 we performed 25000 sweeps at each of R Q = 0,2,5, 10,20, 
with the data divided into 20 bins in each case. At (3 = 21 we performed 20000 sweeps at 
each of R D = 0, 3, 5, 7, 10, with 20 bins in each case. 

The quantities we use in our comparison are, firstly, the masses of the lightest glueballs: 
the ++ , , 2 ++ and, where available, the 2~ + (which should be degenerate with the 2 ++ in 
large volumes and for small lattice spacings). We also use the mass of the lightest flux loop 
that closes through a periodic boundary. This provides us with our estimate of the string 
tension, a, since the mass of this loop is a 2 aL, up to 0(1/ L) finite size corrections, where L is 
the minimal length of the loop (in lattice units). In addition to these masses we also calculate 
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expectation values of the simplest closed loops made out of our 'blocked' links. At a 'blocking' 
level of unity, B t — 1, we have the simple plaquette. More generally these 'superplaquettes' 
consist of a square that is length 2 Bl ~ l in lattice units. The simple plaquette is dominated by 
ultraviolet fluctuations and is of relatively little physical interest. At higher smearing levels, 
the expectation value is dominated by fluctuations closer to physical length scales and how 
the accuracy of these is affected by overrelaxation is a more interesting question. 

We compare the errors on these quantities in the different runs characterised by different 
values of R Q . The reference run is the one with R Q = 0, i.e. pure heat bath. If the R a = 
error is changed by a factor of, say, 7 in the mixed run then the latter is as good as a pure 
heatbath run whose length is I/7 2 times that of the mixed run. This of course assumes that 
our bins are large enough to be essentially independent. We have performed a variety of 
checks to convince ourselves that this is the case for the results we present here. For example, 
for the (3 = 15 pure heatbath run we checked that the bins could be made a factor of 10 
smaller and still be negligibly correlated. (With a factor of 20 the independence began to 
break down.) In order to keep our bin sizes sufficiently large so that we could be confident of 
their mutual independence, the number of bins for each calculation could not be made very 
large. Hence there will be substantial fluctuations on our error estimates. For this reason the 
reader should be cautious about drawing conclusions from any one error ratio, and in practice 
we will average the error comparisons over several quantities. 

In Tables 0, |3] and |4] we show the statistical errors for the flux loop and glueball masses 
at the three different values of (3. There are 2 rows of numbers for each mass. The second 
row contains the actual mass estimates. The first row is the error on the effective mass 
extracted from the same correlation function but from one time step earlier. This contains an 
admixture (typically only a few percent) of excited states. (See Section 3.2.3 for a discussion 
of effective masses.) We display both because the individual error estimates contain quite 
large fluctuations which appear to be largely independent and so can be averaged to obtain 
more reliable error ratio estimates. Because we are equally interested in all these physical 
quantities, it makes sense to construct a global average of these error ratios. We attach to this 
average an 'error' obtained by treating the variations of the individual error ratios around the 
global average as though they were statistical fluctuations. This is intended to do no more 
than provide an indication of the significance of the value of the average error ratio. The 
reader can manipulate the numbers in the Tables in other ways if he so prefers. 

Consider first (3 = 11.0 (Table |2|) and the ratio of errors in the run with overrelaxation 
to the errors in the pure heatbath run. Taking the ratio of corresponding errors in the two 
columns we obtain an overall average error ratio of 1.02(6). So in this case there is no 
improvement in incorporating overrelaxation. At (3 = 15.0 (Table |3|) we obtain average ratios 
1.27(12), 0.96(7), 1.09(9), 1.09(10) for R Q = 2,5, 10,20 respectively. Again there is no sign of 
a significant improvement for any overrelaxation mix. The global error ratio average, 1.10(5), 
confirms this. 

We turn now to our calculation at the weakest coupling, (3 = 21.0 (Table |J). We obtain 
average ratios of 0.91(8), 0.81(6), 0.88(7), 0.94(7) for R D = 3, 5, 7, 10 respectively. We observe 
a clear reduction in the errors of the runs with overrelaxation: the global error ratio average 
is 0.885(33). Although we cannot be certain which mix is best, there is evidence that a ratio 
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of around 5:1 to 7:1 is as good as any at this value of (3 and that this leads to an error ratio 
of around 0.84. To this improvement we should add the fact that an overrelaxed sweep, for 
SU(3) in D = 2 + 1, takes about 77% of the time for a heat bath sweep. Thus the gain in 
using over-relaxation is about 40% in the update time. Although the gain in the total time 
will be reduced by the inclusion of measurements (typically tuned to be about half of the total 
time) there is no doubt that this is a worthwhile gain. 

We turn now to the case of SU(2) where our tests are much more limited. The lightest 
masses here are of the flux loop and of the + and 2 + glueballs. We performed comparisons at 
(3 — 6.0 and (3 = 9.0. Using the calculated string tension to set the physical scale, these values 
of (3 correspond to (3 ~ 2.4 and 2.55, respectively, in the D = 4 SU(2) theory. At (3 = 6.0 
we see no sign of any benefit from overrelaxation; albeit in a calculation of limited statistics. 
Our calculations at (3 = 9.0 (Table ^) are with R Q = 0, 5, 9, 49, 249, oo and have much better 
statistics: 25000 sweeps, split into 25 bins, for each value of R Q . The lattice is 12 2 24 which is 
of small, but not very small, physical size: there are certainly some finite size effects involving 
the 2 + . We see from Table [5] that there is a significant benefit to using overrelaxation. We do 
not show our results for R Q = 249 and R Q = oo (all sweeps overrelaxed) which, while amusing 
for various reasons, are not really relevant to this study. We obtain average error ratios of 
0.76(8), 0.71(9), 0.80(11) for R Q = 5,9,49 respectively. The global error ratio average is 
0.76(5). For SU(2) overrelaxation is a simple operation and is much faster than the heatbath; 
a run with R Q = 5 or 9 is about twice as fast as a run with only heatbath sweeps. Thus the 
overall saving is a factor of 2 x 1/0. 76 2 ~ 4. This is a large reduction. Again, the inclusion of 
measurements will reduce the gain somewhat. 

We turn now to the smeared superplaquettes. The average error ratios, for the runs 
described above, are summarised in Table || We note that, not surprisingly, the simple 
plaquette acquires a larger error if we include overrelaxation. (These global averages mask 
the fact that for small R Q the error is often reduced, which is surprising.) On the other hand 
we observe that the errors on large superplaquettes are reduced and that, in contrast to what 
we saw for masses, this effect is present at smaller couplings. This is similar to what has been 
found in 4 dimensions |29] for large and blocked Wilson loops. 

In conclusion, we have seen that for sufficiently small couplings - equivalent to D = 4 
values of (3 ~ 6.15 for SU(3) and (3 ~ 2.55 for SU(2) - there is a substantial increase in 
efficiency through mixing heatbath and overrelaxation sweeps. The cpu saving is about 40% 
for SU(3) and about 75% for SU(2). The difference is largely due to the fact that SU(2) 
overrelaxation is a very simple and fast operation. In D = 4 the operation of calculating 
'staples' is a little lengthier and so this effect will be somewhat weaker there. 



Appendix B : Testing the efficiency of the operators. 

In constructing a 'good' basis of operators for our various mass calculations, the use of 
spatial blocking is crucial. The general motivation is that if one wants a good overlap onto the 
lightest physical states then one needs to employ (combinations of) large smooth operators. 
There are obviously many possible variants on the particular recipe we have used in this paper 
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(which is the one that has been used successfully in earlier D = 3 + 1 mass calculations). In 
the first part of this Appendix we will consider some variation in the blocking procedure and 
we shall see that our choice is indeed an efficient one. 

The range of operators we have used has been limited not only in the type of blocking 
employed but also in the variety of ways we put the blocked loops together to form colour 
singlet operators. In practice we limited ourselves to simple closed loops. However once we 
go beyond SU(2) there is a whole new class of operators that we can construct, and which 
take advantage of the fact that one can tie together N c indices with a totally anti-symmetric 
tensor. We refer to this, for obvious reasons, as a baryonic vertex. Such operators have not 
been used in previous lattice glueball calculations as far as we are aware. Our attention was 



drawn to them by their possible role in splitting the C = ± sectors, as pointed out in [11 
For this reason we have carried out a small calculation to check whether they encode some 
interesting new information. This is described in the second part of this Appendix. 

Variations on the blocking procedure 

In our construction of 'blocked' link matrices the most obvious parameter is the weighting 
of the direct path as compared to the staple-like paths. The choice we made was to take an 
equal weighting for all the paths. So, for say the x-direction, we would take xx+yxxy^ +y^xxy 
using the notation in Section 3.2.1. In this section we shall perform some calculations using 
a variable weighting 7^. That is to say, we use a blocking 

r'J = i d U x U x + U y U x U x Ul + UlU x U x U y (54) 

where we have suppressed some obvious arguments etc. The blocked link, U B , is then obtained 
by projecting U B back into the group. We then see which value of 7^ is most efficient in the 
sense of producing the best operators. 

Before doing so we briefly comment on the projection back into the group and the resulting 
gauge transformation properties of the U B . We begin by noting that if we perform a local 
gauge transformation on the fields then U B — > g n U B g n i where g n is the gauge transformation 
at site n and the paths making up U B start at the site n and end at the site n'. If the group 
is SU(2) we obtain U B by dividing U B by det{U B }. Since the matrices g are unitary we have 
detg n = 1 and so det{U B } is gauge invariant. Thus in the case of SU(2) U B has the gauge 
transformation properties of a product of links from n to n! and we can form colour singlet 
operators out of closed loops in the usual way. For SU(iV c 7^ 2) the situation is different: U B is 
not proportional to an SU(iV c ) matrix and if we want U B to be in the group we need to define it 
some other way. The method we use is to define U B as equal to the value of the SU(iV c ) matrix 
U that maximises Tr{U B W}. It is easy to see that for SU(2) this reduces to the method we 
use there. It is also trivial to see (using the cyclic property of the trace) that if U B — > g n U B g n , 
then, just as in SU(2), U B — > g n U B g n > and we can form colour singlet closed loops in the usual 
way. However, in practice we maximise the trace by a simple iterative procedure which we stop 
before complete convergence in order to save computer time. This procedure requires, as its 
starting point, some first guess, U s , for the blocked matrix. In practice we construct U s from 
U B in such a way that it does not transform as U s — > g n U s gni under a gauge transformation 
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g. This means that when we stop the algorithm prior to complete convergence, the resulting 
U B only transforms approximately as U B — > g n U B g\ l ,. In principle this does not matter; 
averaging over all field configurations in the Monte Carlo will lead to a cancellation of the 
non-gauge invariant pieces in the correlation functions. However, again in practice, this means 
we generate extra noise and this will increase our statistical errors - something to be avoided 
if possible. 

We see from the discussion in the previous paragraph that there is more to 'blocking' 
than choosing a sum of paths and a relative weighting. One can ask if projecting back to the 
group produces better operators than not doing so (and perhaps using some other form of 
normalisation). Studies in D=3+l of several alternative strategies in SU(2) pi and SU(3) 



32] suggested that this was more-or-less so. Some tests then showed that approximating the 
maximisation of the trace by one or two iterations did not significantly worsen the operators 
or increase the errors. However there has been no demonstration that this continues to be the 
case as we increase the size of the group or that all this continues to hold in D=2+l. These 
are studies that still need to be carried out. 

We return now to our study of how the operators vary with the choice of 7^ in eqn([54]). 
Our calculations are in SU(2) and are performed on a 16 3 lattice at (3 = 7.5 in a run consisting; 
of 10000 sweeps. On these configurations we performed separate mass calculations, using our 
usual basis of operators, for the 5 different blocking schemes that used 7^ = 0.25, 0.50, 1.0, 2.0 
and 4.0. 

What we want to know is how efficient are the different schemes; and in particular how 
efficient is our usual choice, 7^ = 1. We shall confine ourselves to the lightest states in 
each J PC channel. In that case our usual variational criterion, as discussed in Section 3.2.3, 
provides us with a simple criterion for comparing operators of the same quantum numbers: 
one calculates the effective mass at t — a and the 'best' operator is the one that gives the 
smallest value of this effective mass. So what we have done here is to find the best operator 
in our basis for each type of blocking. The best form of blocking will then be the one that 
produces the minimal value of the effective mass. 

In Table [7| we present the value of am e ff(t = a) for the best operator for each of the 5 
kinds of blocking we consider, and for the various J PC quantum numbers. We observe that 
7d G [1, 2] seems to work best overall, although 7^ = 0.5 is virtually just as good if we ignore 
the 1 states. Note that since the different calculations are performed on exactly the same 
sequence of field configurations, the errors will be highly correlated. 

We have therefore seen that with respect to variations in this particular parameter, our 
choice of 7d = 1 is about as good as any. Of course one can vary the algorithm in many other 
ways; for example by including other paths than just the direct path and the 'staples'. A 
systematic study would be useful. 

Operators with 'baryonic' vertices 

We will consider the specific case of SU(3). Suppose we have three curves C%, C 2 , C 3 each 
of which starts at some point n and finish at some point n! . Let us denote by f/ 1 ,?/ 2 ,?/ 3 
the corresponding path ordered products of (blocked) link matrices along these three curves, 
running from n to n' . We can form singlet operators out of pairs of these in the usual way, 
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(55) 



•y - eijkU iif U j:jl U kk ,ei'fk> 



(56) 



where we have exposed the matrix indices and is the usual totally anti-symmetric ten- 
sor. This extends to iV c > 3 in the obvious way; we have U 1 , . . . , U Nc paths joined by the 
appropriate iV c -component e tensor. 

Since the operators in eqn(|55D and eqn(^) have the same quantum numbers they will 
have non-zero overlaps and there is no d priori reason to think that we have lost anything 
by excluding the latter. However it might be that they constitute more efficient operators for 
some states and if that is the case for one of the heavier states, where in practice we cannot 
calculate correlators beyond small t, it might be that, in using them, we will expose a state 
that we have not been able to see with operators of the type in eqn ([55|) . 

In this Appendix we will describe a small exploratory calculation designed to see if includ- 
ing such operators might make a serious difference to our calculations. We shall consider an 
operator, tpy, of the form in eqn(|56"D with 



suppressing obvious arguments and indices. This is a rectangle with a central link crossing 
the rectangle. The path ordering is out from the same vertex for all three curves. Under C 
such an operator reverses all three arrows on the curves; which for this particular operator is 
equivalent to a rotation of tc. So one can easily see that from this operator (and the one we 
obtain by x «-> y) we can obtain ++ , 2 ++ and 1 quantum numbers. 

We have performed a calculation on a 16 3 lattice at (3 = 15 with this operator. In Table |S| 
we list the effective ++ ,2 ++ and 1 masses obtained at t — a using a basis that includes 
the best two blocking levels of this operator. We compare it with what we obtain (on the 
same set of field configurations) if we do not include 0y. We observe that there seems to be 
nothing new in these channels when we include 0y; at least not as far as the ground state and 
the first few excitations are concerned. There is a slight improvement in some of the overlaps, 
as indicated by a decrease in m e ff(t = a), but one would expect that just from an increase in 
the size of the basis. 

This, albeit minimal, study leads us to believe that the inclusion of operators incorporating 
baryonic vertices will not alter our conclusions in any significant way. However such operators 
can be convenient in providing a simple means for constructing J = 1 operators; and they 
may well be important in investigating some physics, e.g. the C — ± splittings, lower order 
in l/N c corrections, . . ., so a more detailed investigation would be useful. 



U 1 

u 2 
u 3 



U X UyUl 
Uy 

UlU y U x 



(57) 
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Appendix C : Testing the benefits of mean-field improve- 
ment 



In this Appendix we shall show that the mean- field/tadpole improved inverse coupling ED 



h = 1 x (jfTriU,)) m 



provides a much better expansion parameter than (3 in our D=2+l calculations. This both 



complements the available D=3+l evidence |50J and provides us with a more accurate way to 
determine lim^o \J° 1 9 1 - 

Our strategy will be to compare directly various extrapolations to a = using either 1/(3 
or as expansion parameters. We shall perform these comparisons using our calculations 
in the SU(2) and SU(3) theories, since these cover large ranges in (3. Having found which 
extrapolation works best, we shall take that information over to the SU(4) and SU(5) theories 
where the range of our calculations is much more limited and where the use of a good expansion 
parameter pays significant dividends. 

We have focussed upon the string tension because in practice this is the quantity that we 
calculate most accurately on the lattice. Since lim a ^ P — 2N c /ag 2 , we know that 

lim (3a Jo = 2N C lim ^ (59) 

(3-+oo a-+0 g z 

The approach to the continuum limit will involve higher order corrections that are inverse 
powers of (3 and hence vanish as powers of a. We thus expect that the approach to the 
continuum limit will be under much better control than in 4 dimensions, where the analogous 
quantity that one would be calculating is lim^o J° I ^mom and where the corrections would be 
inverse powers of log a. This will allow us to make a much more explicit and direct comparison 
than is possible in D=3+l. 
Since we expect 

(3aJa = c + j + j 2 + . . . (60) 

for large (3, it is useful to plot the values of (3aJa against 1/(3] we do this in Fig.|| for SU{2) 
and in Fig|4| for SU(3). For large enough (3 the first two terms in eqnflSTf) will dominate and so 
the values should fall on a straight line as we approach the continuum limit. This we observe 
to be the case. For orientation we also show in Fig|| and Fig.|] the strong coupling predictions 
for the string tension up to 0((3): 

aV = -log(^) + 0(/3 2 ) (61) 

for SU(2), and 

«V = -log(^)-^ + 0(/3 2 ) (62) 

in the case of SU(3). (The extra 0((3) term in the case of SU(3) arises because in that case a 
product of two plaquettes can be used just as well as a single plaquette in tiling the minimal 
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surface spanning the Wilson loop.) We see that our calculated values of the string tension 
extend well into the strong-coupling regime. In this region an expansion in 1//3, such as in 
eqn(|60|) should no longer be valid. 

The maximum range, more-or-less, over which we can perform linear fits with acceptable 
X 2 , turns out to be (3 > 4.5 in the case of 377(2) and (3 > 15.0 for the case of S77(3). These 
fits are given by 

pa^a = 1.324(12) + L2Q ( 11 ) : SU(2) (63) 

P 

and 

Pay/a = 3.275(24) + 8 ' 35 ( 61 ) : SU(3) (64) 

P 

(Note that the errors on the intercept and slope are anti-correlated.) 

While such a linear extrapolation is a perfectly acceptable procedure for extracting the 
continuum value of the string tension, it must suffer from some sytematic bias due to the 
neglect of higher order terms. These, it is clear from the figures, are certainly not negligible 
at intermediate values of (3. If we include 0(1/ f3 2 ) terms in our fits, we naturally find larger 
acceptable ranges for the fits: (3 > 3.0 for 577(2) and (3 > 6.5 for £77(3). The fits are 

(3aV<J = 1-337(23) + + HM ; S U(2) (65) 

(3 (3 l 



pa^/a = 3.367(50) + - \ ' ' + q2 ■ SU(3) (66) 



and 

4.1(1.7) 46.5(11.0) 

T~ + J 2 

We observe that in both cases the inclusion of the extra 0(1/ P 2 ) term has increased the value 
of the continuum limit by an amount that, while small in absolute units, is uncomfortably 
large when compared to the claimed errors, especially so in the case of SU(3). Moreover, in the 
case of SU(3) the coefficient of the 1/P 2 correction is so large that the value of this correction 
is comparable to that of the 1/(3 correction over much of our range. Under such circumstances 
one cannot motivate the neglect of the next, 0(1/P 3 ), correction. However it is clear, from 
the large errors in eqnQ5"6|), that our SU(3) data will not be able to resolve these higher order 
terms with any useful accuracy. Moreover there is also the danger that the 0(1/ P 2 ) correction 
is being overly biased by the values of a^/a in the transition region between weak and strong 
coupling, where the very validity of an expansion in 1/P is breaking down. This leaves us 
with an intrinsic systematic error on the SU(3) continuum limit that may well be larger than 
the quoted statistical error. 

The lattice corrections in eqn(^) are precisely what the use of a better coupling should 
improve - by reducing their coefficients. How well does that work here? If we use eqn(|59|) to 
define Pi we can plot Pidyja against 1/ Pi as in Fig.|| and Fig.|. It is immediately apparent 
from a comparison with Fig.|3| and Fig.[|. that in terms of the 'improved' coupling the higher 
order lattice corrections are dramatically reduced. More quantitatively, if we perform fits as 
before but with Pi replacing P, we obtain the following results. In the case of SU(2) we obtain 
good fits with just the leading 0(1/ Pi) correction for the much larger range (3 > 3.0 while 
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for SU(3) excellent fits are possible for j3 > 8.175 (and reasonable ones all the way down to 
(3 = 6.0). Moreover these fits 



Pidy/a = 1.341(7) 



0.421(50) 



: SU{2) 



(67) 



2.43(22) 

fray/a = 3.318(12) : SU(3) (68) 




display much smaller corrections to the leading asymptotic terms than was the case in eqns(|5"3"D 
and (|64|). Since these fits are so good, there is no real reason to include higher order corrections. 
However if we do so then we obtain 



/W- 1.336(9) -^>-™ :SU(2) 

Pi Pi 



faa^a = 3.323(28) - ^IM. + ^1^1 ; ^(3) (70) 




We observe that the coefficients of the higher order terms are small: so there is no reason 
to worry about the next correction. A second and related observation is that the asymptotic 
values are little changed with the inclusion of the 0(1/ Pi 2 ) correction - in contrast to what 
happened when we used f3 as our expansion parameter. Indeed even the coefficients of the 
terms are insensitive to the inclusion of a higher order term. All this represents a substantial 
improvement in the perturbative control of the continuum limit. 

The fact that we can extrapolate fiia^/a with fits involving just two parameters, means 
that we do not need to perform calculations at more than four values of f3 in the case of SU(4) 
and SU(5). This represents a substantial saving in computational effort. 

From fits such as the above we can extract the continuum mass ratios shown in eqn(^). We 
remark that it is both because we are in 3 dimensions, where the bare coupling decreases lin- 
early with the scale a rather than just logarithmically, and because of the extent and accuracy 
of our lattice calculations, that it is possible to perform reasonably accurate extrapolations to 
the continuum limit even with the 'bad' lattice bare coupling. This has enabled us to quantify, 
in a way that is not yet possible in 4 dimensions, how much the mean-field improved coupling 
actually improves the approach to the continuum limit of the lattice spacing, a. 

Appendix D : Calculations with an asymmetric lattice 
action 

In this Appendix we present an SU(2) calculation of the mass spectrum on lattices with 
timelike and spacelike lattice spacings related by a t ~ a s /4. As discussed in Section 5.1 the 
primary purpose of this study is to check explicitly that our criteria for which effective masses 
adequately reflect the actual masses, are in fact accurate. A second reason is that this, being 
a calculation with a different action to the one we have used so far, will provide us with 
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some test of universality. We shall first discuss some of the features that are peculiar to such 
calculations. We then present our results. 

Preliminaries 

To allow different spatial and temporal lattice spacings we use the action in eqn@. What 
is the relation between our choice of (3 S , {3% and the lattice asymmetry? Suppose we are aiming 
for a particular ratio 

r = X (71) 



In the limit a s , a t — > we have 



{l--ReTrU Pa} -+ at—TrF 
{1 - ±-ReTrU Pt } -> a y t ^-TrF 2 

(3S -> 4 / d 2 xdt-TrF 2 (72) 

where F 2 is the continuum field strength squared. Since the integration measure gives a factor 
a 2 a t when discretised, we see that the choice 

(3s = r(3 ; & = /3 = ^ (73) 

r a s g 2 

is what is needed, at least at the classical level, to achieve the asymmetry, r, defined in 
eqnfl7ip. 



In practice there will be quantum corrections to these classical relations. Three related 
questions immediately arise. 

(a) What do we need to know about a s and ct£? 

(b) In a given simulation, how can we calculate a s and at directly? 

(c) Can we easily 'improve' upon the relations in eqn(|73|)? 

Before considering each of these questions in turn, we need to remark that the classical 
relations in eqn([73]) should remain a roughly reliable guide in the full theory. This is because 
our theory is super-renormalisable and this is in contrast to the situation in 4 dimensions. 
Nonetheless we do expect significant corrections, as we saw when considering the (3 dependence 
of (3ayja. We saw there (see Appendix C) that the corrections to the classical relation (3 = 
2N c /ag 2 are quite large and can be drastically decreased by the use of a mean-field improved 
coupling. 

In the rest of this Appendix we shall assume that the asymmetric lattice has been chosen 
so that a t <C a s . This means that we shall systematically ignore any 0(af) corrections as 
compared to to ones that are 0(a 2 ). 

(a) How well do we need to know a s and a{l 

The first thing we need to establish is whether we actually need to know a s and at any 
more accurately than we already know them through using eqn(|73|). 
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If we just wanted to calculate some lattice mass ratios, mj/m , then these could be obtained 
without knowing the lattice spacing at all: our usual procedure would give us estimates of 
a t mi and the lattice spacing then cancels in the ratio. However if we want to extrapolate to 
the continuum limit then the leading correction will be 0{a 2 ), assuming at <C a s as will be the 
case here, and so we need to know a s m for some mass m in order to provide the correction term 
in the analogue of eqn(pE9|). This can, however, be finessed by using eqn(|50D instead; although 
1//3 2 ~ a 2 g 4 /4N 2 + 0(a 3 ), the correction has a small enough coefficient that it should not 
significantly degrade the accuracy of our extrapolations. 

Of course, we also need to control finite volume effects in an efficient way - that is to say, 
more efficient than doing a detailed finite volume study at each coupling. So we need to be 
able to compare the lattice size, L = a s L s , at different couplings. For this we need to know 
a s to a reasonable approximation; great accuracy is not needed because we usually include a 
margin of safety in our choice of the volume. 

There is however at least one place where we do need accurate values for the lattice 
spacings. This is in our calculation of the string tension. Our usual procedure is to calculate 
the mass, a t mp(L), of a flux loop of length L = a s L s that winds around the spatial torus. 
This mass, a t mp(L), can be written, using eqn(|32|), as 



a t mp(L) = a t x La(L) 
= a t L[a(oo) - 



7T 



6L 2 

T 71 at (>7A\ 

= a t a s L s a - — (74) 

bL s a s 

Clearly we need to know r = a t /a s very accurately if we are to be able to calculate a t \J<J with 
the accuracy we are used to. A similar situation arises if we calculate potentials using Wilson 
loops. 

(b) How do we calculate a t /a s directly? 

There are two obvious methods that we can use to calculate r = a t /a s . The first involves 
calculating the energies of states with non-zero momenta. Suppose we have a particle of mass 
m. The allowed momenta are a s p = (2irn x / L s ,2irn y /L s ) and the corresponding energies that 
we obtain from our correlators may be written as atE(p). For small momenta we expect the 
continuum dispersion relation, E 2 = p 2 + m 2 , to be accurately satisfied. (We have explicitly 
seen that this is so on the symmetric lattices that we have used in our main calculations in 
this paper.) So we expect to have 

a 2 E 2 (p) = a 2 (p 2 + m 2 ) 



-y(^) 2 (n 2 x + n 2 y ) + (a t m) 2 (75) 

da Li 



Therefore, from our calculated values of a t E(p) and a t m we can obtain, using eqn([7^), a value 
for a t /a s . 

Our second method is even more direct. Normally we calculate correlators in the t- 
direction. We could instead calculate our correlators in, say, the x-direction. In that case 
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our space would be (y,t) in place of (x,y). As long as a t L t is large, as it will always be, this 
new spatial volume will also be large and we can assume there are no finite volume effects. 
Thus we obtain the same mass in both calculations, up to lattice spacing corrections. That 
is to say, we obtain ajm from our t-correlators and a s m from our ^-correlators. Equating the 
two we obtain a t /a s . 

We can combine the above two methods by calculating glueball correlation functions in 
the x-direction with non-zero momenta in either the y or t directions. Comparing the energies 
of such states gives us another direct estimate of a t /a s . If we include a range of momenta 
one can attempt to tune this ratio so that the dispersion relations in p t and p y coincide. This 
allows us to fix the asymmetry without assuming the continuum dispersion relation. 

In practice, we shall not use this last method, and we shall only consider the lowest two 
momenta in applying our first method. While this reduces the precision with which we can 
estimate the ratio a t /a s , it suffices for our purposes. 

Both of the above methods will suffer from lattice spacing corrections. The continuum 
dispersion relation will only be valid up to corrections of order (a s p) 2 and the eigenstates of 
the transfer matrices defined on the x, y and y, t spatial tori will differ by order a 2 s corrections 
because the latter torus has a lattice spacing a t rather than a s in one of the two directions. 
However as long as we are consistent in the method used to estimate a t /a s we can absorb this 
correction into the correction term used in taking the continuum limit. 

Before turning to some explicit calculations of the above kind it is worth pointing out 
that although the second method described above seems more direct, it is in practice more 
awkward to implement. The reason is that we need to produce blocked link matrices in order 
to have useful operators, and this has to be carefully tailored in the case where one spatial 
lattice spacing is very different from the other. And this is in addition to the fact that using 
two different spatial planes means producing two sets of blocked links. For these reasons our 
calculations using the second method will be on only a subset of our lattices. 

(c) Can we 'improve' upon the estimate of r? 

We have seen in our previous calculations that we get much smaller corrections to the 
limit (3a^/a — > 2N c ^/a/g 2 as a — > 0, if we use a mean-field improved inverse coupling, in 
place of (3. One might hope that a similar approach with an asymmetric lattice action would 
improve our control over the value of r. The straightforward implementation of this idea in 
the context of the action in eqn(|8|) would be to define 'improved' values of r and f3 by 

riPi = r(3{^TrU Ps ) 



This is in the spirit of the approach suggested in JE| although there one effectively replaces 
jj-(TrU Pt ) by unity. We shall calculate 77 below and explicitly check how much of an improve- 
ment one really obtains. 
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The calculation 



We perform calculations at (3 = 4.0, 5.3 and 8.0. In all three cases we choose r = 0.25 in 
the action, as given in eqn(|8]). The lattice sizes are 12 2 60, 16 2 64 and 24 2 96. If the classical 
relations in eqn([F3|) were valid, then the value of a s would be exactly what we obtained at the 
corresponding values of (3 on symmetric lattices. The reader will note that our lattice sizes 
are somewhat larger than would be necessary if this were the case; this is to give us some 
margin in case the quantum corrections to these relations are significant. 

In addition to these calculations we also perform calculations on somewhat smaller lattices, 
8 2 60, 12 2 64 and 16 2 96 respectively. It is on these lattices that we calculate correlators in both 
t and x directions. We shall only calculate the mass of the periodic flux loop, mp(L), on these 
lattices. Smaller lattices are preferable for this purpose because mp(L) will be smaller and 
so we will have more accurate values. (Obviously this is only important for the correlators in 
the x-direction where a s mp will not be small.) 

Our original motivation for performing such calculations was to have a finer resolution 
on the effective mass plot, so as to see whether the typical heavier mass could really be 
extracted from its effective masses between t — a to t — 2a. Since we have a t ~ a s /4, this 
question becomes: have the effective masses in the region Aa t <t< 8a t already reached their 
asymptotic plateau? So without further ado we plot the effective masses from the 24 2 96 lattice 
in Fig.|TC| Note that if we wished to obtain such a small lattice spacing in our usual symmetric 
lattice calculation, we would have had to do it on a 96 3 lattice at (3 ~ 30: a daunting prospect! 
We see from Fig.|TU| that there is reasonably good evidence in all cases that an effective mass 
extracted over the interval Aa t < t < 8at provides an unbiased estimate of the asymptotic 
mass. We infer from this that extracting masses from the range t = a to t = 2a on symmetric 
lattices in the neighbourhood of /3 — 8, where the heavier states are in the noise for t > 2a, 
is in fact justified. Our j3 = 5.3 calculations also support this way of calculating masses 
(although with less precision) thus reassuring us that the estimates we have used in this paper 
are indeed unbiased over the whole range of f3 relevant to our continuum extrapolations. 

Since these lattice actions are different from the symmetric ones that we have used in the 
body of the paper - indeed, as we have already remarked, one may regard them as being close 
to the Hamiltonian limit - it is interesting to extract a continuum mass spectrum from them, 
so testing universality, to some extent. 

In Table |3D] we list the masses we have extracted at our three value of (3. They are all 
in units of a t . In addition to the glueball masses we also list the mass of the flux loop that 
winds around the spatial torus. As discussed above, we need to know the value of a t /a s in 
order to extract a ty Ja (or a s y/a). And we need to know a ty /a if we are to calculate m G /^/a 
for comparison with our previous calculations. We therefore turn to this next. 

In order not to confuse different quantities, we shall continue to use r for the parameter in 
the action, and we shall choose r = 0.25 here. Classically, but only classically, we know that 
r = a t /a s . The 'true' value of a t /a s is the one that we explicitly calculate using the methods 
described earlier in the Appendix: this we shall either label r meas or simply refer to as a t /a s . 
Finally there is 77 as defined by eqn([75|). 

We calculate a t /a s by comparing the flux loop energy as calculated from two momenta, 
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Pi and P2, and then using eqn([75]). We do so using the lowest three momenta which we shall 
refer to as Pi — 0, 1, 2 for convenience. This we shall do on both our larger and smaller lattices 
at each value of (3. Note that the smaller p the more reliable will be the continuum energy- 
momentum dispersion relation. Since the lowest momenta decrease as the lattice size increases, 
the values obtained on the larger lattices will have smaller systematic errors although larger 



statistical errors. The values thus obtained are listed in Table 31. We also show there the 
values we get for a t mp/a s mp; these have only been calculated for the smaller lattices. From 
these results we infer the following values of the asymmetry: 



a t 



0.245(5) (3 = 8.0 

0.237(6) (3 = 5.3 (77) 

0.230(10) /3 = 4.0 



The deviation from the classical value, a t /a s — r — 0.25, is not large. We also note that 
although the 'improved' value, 77, is shifted in the right direction, it overshoots so that it 
is no closer to the 'true' value than is r. This is even more so if one replaces the time-like 
plaquettes by unity in eqn([76]). Thus it seems that the most naive mean-field 'improvement' 
is not an improvement here. 



Using in eqn (|74|) the values of a t /a s in eqn(|77|) and the values of a t m P in Table [30], we 
obtain 

f 0.04408(56) (3 = 8.0 

a^a = I 0.06806(102) (3 = 5.3 (78) 

[ 0.0919(21) (3 = 4.0 

The uncertainty in a t /a s has roughly doubled the error on a ty J<j\ thus it is no longer the most 
accurately calculated quantity (as it was in case of the symmetric action) and if we were to 
calculate mass ratios from scratch we might prefer to use the scalar glueball mass as our basic 
scale. We now extrapolate these values (multiplied by (3) to the continuum limit, using a 
1/(3 correction just as in eqn(P^|). The fit has a very good confidence level and gives us the 
continuum value of Aatyjcr / a s g 2 . In the continuum limit a t /a s = r = 0.25 and so we finally 
obtain 

^ = 0.3375(130). (79) 
T 



This is certainly consistent with the value of 0.3353(18) in eqn(|38D which was obtained with 
r = 1, but the error is much larger. In large part this is just because the present calculation 
is a much smaller one. But in some part it is due to our uncertainty in the value of at/a s . 
Without this uncertainty our error in eqn(|7^) would have been smaller by about a factor 
of 1.5. We have performed other continuum extrapolations as well. If we use the classical 
value a t /a s = r = 0.25 in our calculations, then we again obtain a good fit, this time with a 
continuum value of 

^ = 0.3288(50). (80) 
T 

This is consistent with our previous values, as it should be because the corrections to r can be 
absorbed into the 1/(3 correction. Indeed we find that the coefficient of the 1/(3 term is larger 
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in the latter case: ~ 0.22(3) versus ~ 0.12(8) when we use the values in eqn([77|). We have 
performed other extrapolations as well: the errors vary but the values are consistent with each 
other. It is worth remarking that the fit leading to eqn(|79|) has a much smaller correction 
than one finds in the symmetric case when using /?; indeed it is about the same size as one 
obtains using the improved coupling, /?/. This suggests that lattice corrections are smaller on 
very asymmetric lattices, and perhaps explains why we did not gain anything from using the 
mean-field recipe. 

We can now take the values of a t m G in Table [Jt], the values of a ty /a in eqn([78|), form 
ratios, and extrapolate to the continuum limit using eqnflSOl). We obtain the continuum mass 
ratios shown in Table |32|. We also show there the symmetric lattice values that have been 
obtained elsewhere in this paper. We observe that they are consistent within errors. We note 
that if we extrapolate with an 0(a 2 a) correction as in eqn([|5p, we obtain almost identical 
results. This is also the case if we use string tensions calculated using a t /a s = 0.25, except 
that the fits tend to be significantly worse - as one would expect if this involved an error that 
was really 0{lj 0) rather than 0(1/ (3 2 ). 

We observe in Table 32 that the errors on the lighter masses are larger in the asymmetric 
case. This is no surprise since our symmetric calculations are very much larger. What is 
striking is that for the heaviest masses, such as the 1 ±+ , the asymmetric errors are actually 
smaller. This displays the power of such calculations for determining the masses of heavier 
states. 

Two final asides on the spectrum. For reasons we do not entirely understand, we seem to 
have no problem in obtaining a set of excited 2 ±+ masses that continue well to a = and, 
indeed, are consistent with what one finds for higher groups - see Table This is in contrast 
to the case of a symmetric lattice action. We also note in Table |3(] that the mass of the 
1 h is larger than that of the 1 h *. Our ordering of the states is determined by the effective 
mass at t — a, which in this case seems to be unreliable. This is undoubtedly related to the 
peculiarities in the J = 1 sector that we have previously noted during our symmetric action 
calculations. Obviously it is hard to argue with using the lower mass for the ground state 
1 h , and if we do so we obtain the value in square brackets in Table It is amusing that 
this value fits better with parity doubling and with the values obtained for the higher groups. 
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Appendix E : The U(l) mass spectrum 

In this Appendix we calculate the mass spectrum in the D=2+l U(l) theory. U(l) is as 
far as one can get from U(oo), so this should provide us with a useful contrast to the SU(iV c ) 
mass spectra which we have calculated in this paper. 

Calculating the U(l) spectrum might seem pointless; the continuum limit should be a 
theory of free, non-interacting photons. While this would certainly provide a contrast to our 
SU(iV c ) spectra, it would hardly be very illuminating. 

Although the continuum limit is indeed trivial, in the sense that there will be no bound 
state whose mass is finite in units of the mass scale g 2 , there is nonetheless interesting dynamics 
at finite values of the lattice spacing, a. This arises from the presence of magnetic monopoles 
in the theory. (To be more precise, in D=3+l these would be magnetic monopoles. Here 
they are pointlike instantons whose fields are identical to the spatial fields of a static Dirac 
monopole. Hence we shall follow the usual custom and refer to them as magnetic monopoles, 
even in D=2+l.) These monopoles clearly have an action that is ~ Cul&g 2 where Cm is 
a constant that depends on the particular lattice action being used. Thus a monopole has 
a weighting ~ exp{— cu/ag 2 } and hence the average distance between monopoles will be 
dm ~ aexp{cM/3ag 2 }. (Up to weakly varying factors that come from integrating small 
fluctuations around the monopoles.) This provides a scale for the theory that is different from 
g 2 . This scale is interesting because the monopoles change the physics in a qualitative fashion. 
As is well known [P| they produce a linear confining force between external static charges. 
One could also expect them to produce a nontrivial mass spectrum. At the very least, there 
will be a massive 'photon'. Of course, in the continuum limit d^ — + oo in units of 1/g 2 and 
so on the latter scale the monopoles disappear to r = oo as we approach the continuum limit. 

Actually, the above description represents an over-simplification. There is not just one new 
scale introduced by the monopoles. There is also a scale associated with the string tension, 
l a = l/y/a, and a scale associated with the screening mass, l s = l/m s . These scales mutually 
diverge in the continuum limit: l s ~ I 2 ~ d^ U P to constants and powers of (3. The origins 
of this peculiar situation lie in the fact that the monopoles are singular objects. 

For these reasons we will not try to calculate a 'continuum' mass spectrum. Rather we shall 
calculate the spectrum for lattice spacings that are small compared to the obvious dynamical 
length scale, a <C 1/y/cr. 

The results of our calculation are presented in Table How do they compare to the 
SU(iV c ) spectra listed in Tables and An immediate difference with all of the SU(iV c ) 
spectra is that the ++ is no longer the lightest particle; the is about half the mass. We 
also note that while the ratio m ++/ v /cr was increasing as N c decreased, the U(l) value is 
about as much below the SU(oo) value as the SU(2) ratio is above. Apart from these striking 
differences, the rest of the spectrum seems quite similar (albeit within the large errors). This 
is particularly so if we compare to the ++ mass rather than to ^Jo. For example, the 2 ++ to 
++ mass ratio is close to its SU(iV c ) value. We note also that we have approximate parity 
doubling for J 7^ 0, thus confirming that a is indeed small enough for continuum rotational 
symmetry to have been restored on hadronic length scales. 
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cools 


uq = 4 


n G = 5 


TlQ = 6 


n G = 8 





0.2010 


0.2085 


0.2072 


0.2088 


1 


0.1272 


0.1199 


0.1067 


0.0874 


2 


0.0813 


0.0662 


0.0494 


0.0266 


5 


0.0596 


0.0417 


0.0254 


0.0073 


10 


0.0516 


0.0328 


0.0168 


0.0028 


15 


0.0483 


0.0291 


0.0132 


0.0015 


20 


0.0463 


0.0269 


0.0111 


0.0010 



Table 1: Average action per plaquette when a thermalised SU(5) field is cooled using n G 
SU(2) subgroups. 



state 


R o = 


Ro = 5 


flux loop 


0.0054 


0.0050 




0.024 


0.030 


0++ 


0.0098 


0.0097 




0.056 


0.051 


0~ 


0.019 


0.024 




0.22 


0.20 


2++ 


0.036 


0.031 




0.29 


0.32 



Table 2: Errors on SU(3) masses on a 12 2 16 lattice at f3 — 11; R a is the number of over-relaxed 
sweeps for every heatbath sweep. 



state 


R o = 


Ro = 2 


Ro = 5 


R = 10 


R = 20 


flux loop 


0.0036 


0.0053 


0.0033 


0.0044 


0.0043 




0.0085 


0.0123 


0.0078 


0.0074 


0.0094 


0++ 


0.0086 


0.0132 


0.0075 


0.0139 


0.0085 




0.034 


0.040 


0.037 


0.038 


0.033 


0~ 


0.017 


0.026 


0.020 


0.018 


0.022 




0.081 


0.113 


0.086 


0.082 


0.125 


2++ 


0.025 
0.090 


0.015 
0.092 


0.015 
0.092 


0.024 
0.079 


0.021 
0.070 



Table 3: Errors on SU(3) masses on a 12 2 16 lattice at (5 — 15. R Q is the number of over- 
relaxed sweeps for every heatbath sweep. 
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state 


^0 = 


R Q — 3 


#0 = 5 


Ro = 7 


R = 10 


flux loop 


0.0053 


0.0044 


0.0040 


0.0038 


0.0030 




0.0089 


0.0075 


0.0070 


0.0057 


0.0069 


0++ 


0.0085 


0.0072 


0.0058 


0.0056 


0.0095 




0.0170 


0.0172 


0.0138 


0.0132 


0.0142 


0~ 


0.0123 


0.0072 


0.0090 


0.0150 


0.0103 




0.030 


0.035 


0.035 


0.030 


0.034 


2++ 


0.0148 


0.0106 


0.0103 


0.0119 


0.0120 




0.044 


0.029 


0.032 


0.041 


0.045 


2-+ 


0.0094 


0.0132 


0.0103 


0.0102 


0.0094 




0.043 


0.044 


0.026 


0.044 


0.054 



Table 4: Errors on SU(3) masses on a 24 3 lattice at /3 — 21. R is the number of over-relaxed 
sweeps for every heatbath sweep. 



state 


R o = 


Ro = 5 


Ro = 9 


Ro = 49 


flux loop 


0.004 


0.003 


0.002 


0.002 




0.006 


0.004 


0.003 


0.004 


0+ 


0.007 


0.004 


0.006 


0.007 




0.015 


0.014 


0.009 


0.009 


2+ 


0.011 


0.007 


0.010 


0.010 




0.019 


0.019 


0.017 


0.021 



Table 5: Errors on SU(2) masses on a 12 2 24 lattice at f3 — 9. R a is the number of over-relaxed 
sweeps for every heatbath sweep. 





SU(3);f3 = 21 


SU(3);l3= 15 


SU(3);f3= 11 


SU(2);(3 = 9 


1 


1.28(15) 


1.00(13) 


1.29 


1.01(14) 


2 


0.92(5) 


0.63(6) 


0.88 


1.02(5) 


3 


0.64(5) 


0.63(4) 


0.81 


0.78(3) 


4 


0.75(4) 


0.69(4) 


0.79 


0.71(6) 


5 


1.01(10) 









Table 6: Ratio of errors with and without over-relaxation; for 'plaquettes' at blocking level 
Bl 
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min^, {am 


eff(t = a)} 






state 


Id = 0.25 


Id = 0.50 


Id = i.o 


Id = 2.0 


Id = 4.0 


flux loop 


0.622(11) 


0.612(10) 


0.603(10) 


0.601(9) 


0.630(8) 


0++ 


0.979(15) 


0.967(13) 


0.963(12) 


0.969(11) 


1.008(11) 


2++ 


1.53(4) 


1.55(4) 


1.56(3) 


1.56(5) 


1.63(3) 


2~+ 


1.71(4) 


1.61(3) 


1.58(3) 


1.58(3) 


1.63(3) 


o-+ 


2.51(6) 


2.13(4) 


1.99(3) 


2.01(4) 


2.12(5) 


1++ 


2.80(6) 


2.32(7) 


2.19(5) 


2.20(4) 


2.34(4) 


1-+ 


3.01(19) 


2.43(8) 


2.27(8) 


2.30(6) 


2.48(6) 



Table 7: Effective masses at t — a from the 'best' operators in different blocking schemes, 
described in Appendix B. 





am cff (t = 


; a) 


state 


with <pY 


without 0y 


0++ 


1.08(2) 


1.08(2) 


0++* 


1.67(5) 


1.67(5) 


Q++** 


2.03(4) 


2.05(3) 


Q-| — h*** 


2.26(6) 


2.26(5) 


2++ 


1.78(2) 


1.78(2) 


2 ++* 


2.15(6) 


2.15(6) 


2H — h** 


2.33(7) 


2.34(7) 


2++*** 


2.56(7) 


2.61(7) 


\ — 


2.60(6) 


2.60(6) 


1~ * 


2.62(5) 


2.63(5) 


■y 


2.74(8) 


2.78(9) 





2.92(10) 


3.02(7) 



Table 8: Effective masses at t — a obtained from bases with and without the operator, 
which contains two 'baryonic, vertices. From a 16 3 lattice at f3 — 15 in SU(3). 
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p 


L, V 


°eff 


12 


32,48 


0.66(20) 


9 


24,32 


0.79(26) 




16,24 


0.63(9) 




12,16 


0.47(7) 




8,12 


0.10(2) 




6,8 


0.00(2) 


6 


24,32 


-1.7(2.9) 




16,24 


1.26(70) 




12,16 


0.87(27) 




8,12 


0.41(8) 




6,8 


0.28(5) 



Table 9: Coefficient of effective leading 1/L correction in the flux loop mass: as extracted 
from loops of length L and V using eqn(|3"3"|). 





L 




am(a) 


am(2a) 


am(3a) 


am (4a) 


am(5a) 


14.5 


40 




1 


0.373(1) 
0.375(1) 


0.364(2) 
0.365(2) 


0.364(2) 
0.363(2) 


0.364(3) 
0.361(3) 


0.363(5) 
0.358(4) 


12.0 


48 




1 


0.685(2) 
0.690(2) 


0.668(5) 
0.670(4) 


0.656(9) 
0.671(7) 


0.628(16) 
0.677(13) 


0.632(36) 
0.677(25) 




32 




1 


0.443(2) 
0.447(1) 


0.434(2) 
0.434(2) 


0.429(3) 
0.430(3) 


0.426(4) 
0.425(4) 


0.426(7) 
0.425(6) 


9.0 


32 




1 


0.846(3) 
0.846(2) 


0.826(5) 
0.823(4) 


0.832(10) 
0.821(7) 


0.837(21) 
0.832(17) 






24 




1 


0.617(2) 
0.622(1) 


0.605(4) 
0.609(2) 


0.601(6) 
0.600(5) 


0.590(9) 
0.599(10) 


0.593(17) 
0.597(17) 




16 




1 


0.389(2) 
0.398(2) 


0.381(2) 
0.383(2) 


0.378(4) 
0.378(4) 


0.378(5) 
0.370(5) 


0.379(7) 
0.364(9) 




12 




1 


0.274(2) 
0.291(3) 


0.269(2) 
0.274(6) 


0.268(3) 
0.261(10) 


0.271(3) 
0.268(14) 


0.272(4) 




8 




1 


0.174(1) 
0.215(6) 


0.172(1) 
0.182(10) 


0.172(1) 
0.163(25) 


0.172(2) 


0.173(2) 




6 




1 


0.130(1) 
0.093(16) 


0.129(1) 


0.129(1) 


0.128(2) 


0.128(2) 



Table 10: Effective loop masses for various loop lengths L and from the lowest 2 momenta, 
p = 2im mom / L, in SU(2) at the values of (3 shown. 
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N c 


P 


L 


amp 


C eff 


5 


44 


8 


0.4490(33) 








12 


0.7314(71) 


0.55(8) 






16 


1.0110(56) 


0.75(23) 




33 


8 


0.997(8) 








12 


1.551(25) 


0.53(26) 


4 


28 


8 


0.4257(45) 








12 


0.7139(66) 


0.72(9) 






16 


0.9857(66) 


0.70(23) 


3 


15 


8 


0.4425(27) 








12 


0.7391(46) 


0.72(6) 






16 


1.0103(92) 


0.51(23) 






24 


1.5636(180) 


0.92(44) 



Table 11: Flux loop masses as a function of the loop length, L, for SU(3), SU(4) and SU(5). 
Also shown is c e //, the coefficient of the 1/L correction in eqns (|32] , |33|) . 



p 


L 


a\Ja 


P 


L 


a\Ja 


14.5 


40 


0.09713(20) 


6.0 


16 


0.2538(10) 


12.0 


48 


0.1179(8) 


5.0 


16 


0.3129(20) 


12.0 


32 


0.1179(5) 


4.5 


12 


0.3527(30) 


9.0 


32 


0.1622(4) 


3.75 


8 


0.4487(33) 


9.0 


24 


0.1616(6) 


3.47 


8 


0.4889(56) 


6.56 


24 


0.2297(10) 


3.0 


6 


0.584(16) 


6.0 


32 


0.2529(33) 


2.5 


4 


0.709(11) 


6.0 


24 


0.2562(26) 


2.083 


4 


0.852(50) 



Table 12: SU(2) string tensions as extracted from flux loop masses of length L using eqn(B3). 



p 


L 


ay/a 


P 


L 


ay/a 


34.0 


40 


0.10379(26) 


8.175 


6 


0.5598(31) 


28.0 


32 


0.12753(20) 


7.5 


6 


0.591(23) 


21.0 


24 


0.17479(38) 


7.5 


4 


0.633(3) 


15.0 


24 


0.2570(15) 


7.0 


4 


0.698(5) 


15.0 


16 


0.2553(12) 


6.5 


4 


0.782(26) 


11.0 


12 


0.3748(23) 


6.0 


4 


0.835(41) 



Table 13: SU(3) string tensions as extracted from flux loop masses of length L using eqn(j32l). 
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SU(4) 


SU(5) 


P 


L 




P 


L 


ay/a 


51.0 


32 


0.12859(23) 


82.0 


32 


0.12715(27) 


40.0 


24 


0.16804(30) 


64.0 


24 


0.1664(4) 


28.0 


16 


0.2523(8) 


44.0 


16 


0.2554(7) 


21.0 


12 


0.3597(24) 


33.0 


12 


0.3645(29) 



Table 14: SU(4) and SU(5) string tensions as extracted from flux loop masses of length L 
using eqn(|3"2"D. 



p 


L 


plaq 


P 


L 


plaq 


14.5 


40 


0.929803(3) 


6.0 


16 


0.824744(33) 


12.0 


48 


0.914824(3) 


5.0 


16 


0.786850(20) 


12.0 


32 


0.914823(3) 


4.5 


12 


0.760841(45) 


9.0 


32 


0.885445(5) 


3.75 


8 


0.706986(48) 


9.0 


24 


0.885438(7) 


3.47 


8 


0.680058(59) 


6.56 


24 


0.840548(22) 


3.0 


6 


0.624023(62) 


6.0 


32 


0.824772(10) 


2.5 


4 


0.54737(20) 


6.0 


24 


0.824782(16) 


2.083 


4 


0.47100(13) 



Table 15: Average SU(2) plaquette values. 



p 


L 


plaq 


P 


L 


plaq 


34.0 


40 


0.919680(2) 


8.175 


6 


0.620730(45) 


28.0 


32 


0.901903(2) 


7.5 


6 


0.57810(19) 


21.0 


24 


0.867671(5) 


7.5 


4 


0.57801(8) 


15.0 


24 


0.810773(9) 


7.0 


4 


0.54118(11) 


15.0 


16 


0.810767(11) 


6.5 


4 


0.50060(13) 


11.0 


12 


0.733401(18) 


6.0 


4 


0.45736(19) 



Table 16: Average SU(3) plaquette values. 
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SU(4) 


SU(5) 


P 


L 


plaq 


P 


L 


plaq 


51.0 


32 


0.898791(1) 


82.0 


32 


0.899245(1) 


40.0 


24 


0.869608(3) 


64.0 


24 


0.869510(3) 


28.0 


16 


0.809339(9) 


44.0 


16 


0.805322(6) 


21.0 


12 


0.737628(18) 


33.0 


12 


0.731640(17) 



Table 17: Average SU(4) and SU(5) plaquette values. 



group 


c 


Ci 


CL(%) 


P> 


SU(2) 


1.341(7) 


-0.421(51) 


60 


3.0 


SU(3) 


3.318(12) 


-2.43(22) 


90 


8.175 


SU(4) 


6.065(32) 


-7.74(1.10) 


90 


21.0 


SU(5) 


9.657(54) 


-21.4(2.7) 


70 


33.0 



Table 18: Continuum extrapolations of Pia^/a — > 2N cy /a/g 2 as in eqn(|37D, with confidence 
level of best fit, and range of (3 fitted. 
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state 


am(a) am{2a) am(3a) am(4a) am(5a) 


0++ 


0.541(3) 0.533(5) 0.539(7) 0.523(11) 0.533(14) 


o++* 


0.821(3) 0.798(7) 0.796(14) 0.799(29) 0.795(65) 


Q++** 


1.033(4) 1.009(10) 1.028(28) 1.013(66) 1.18(22) 




0.793(3) 0.779(7) 0.765(10) 0.750(30) 0.71(5) 


* 


1.016(4) 1.002(8) 0.969(26) 1.054(62) 1.09(17) 


Q ** 


1.202(5) 1.193(14) 1.14(4) 1.20(18) 


o-+ 


1.208(5) 1.155(11) 1.146(43) 1.02(11) 


o+- 


1.326(6) 1.256(14) 1.256(59) 1.59(26) 


2++ 


0.913(3) 0.891(7) 0.905(19) 0.911(40) 


2 ++* 


1.103(4) 1.072(9) 1.074(31) 1.00(10) 


2"+ 


0.906(3) 0.879(7) 0.853(14) 0.847(43) 


2 -+* 


1.094(4) 1.076(9) 1.069(32) 1.04(10) 


2- 


1.081(4) 1.054(10) 1.051(28) 1.10(10) 


2 ~* 


1.271(5) 1.225(16) 1.218(46) 1.17(11) 


2+" 


1.075(4) 1.043(9) 1.049(30) 1.10(11) 


2 +-* 


1.276(4) 1.230(17) 1.218(49) 1.36(17) 


1++ 


1.344(4) 1.294(14) 1.284(48) 1.14(18) 


1-+ 


1.355(5) 1.300(10) 1.327(46) 1.33(16) 


1— 


1.255(3) 1.212(9) 1.176(34) 1.18(10) 


1+- 


1.286(5) 1.255(14) 1.231(44) 1.15(19) 



Table 19: Effective masses for the states shown, on a 32 3 lattice at (5 = 82 in SU(5). 
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SU(2) ; P = 6 


state 


L = 6 


L = 8 


L = 12 


L = 16 


L = 24 


L = 32 


0++ 


0.99(2) 


1.12(2) 


1.21(2) 


1.19(2) 


1.20(2) 


1.18(2) 


2++ 


0.92(2) 


1.17(2) 


1.71(5) 


1.80(8) 


1.99(12) 


1.87(10) 


2"+ 


2.20(13) 


2.17(15) 


2.15(21) 


1.81(14) 


2.01(16) 


1.91(11) 


2xl P 


0.62(1) 


0.90(1) 


1.43(1) 


2.00(2) 


3.12(6) 


4.11(15) 



Table 20: How the lightest SU(2) glueball masses depend on the spatial volume. Twice the 
mass of the periodic flux loop is also shown. 



SU(3) ; P = 15 



state 


L = 8 


L = 12 


L = 16 


L = 24 


0++ 
++ 


0.99(4) 
0.98(2) 


1.09(2) 
1.49(4) 


1.15(2) 
1.92(9) 


1.10(2) 
2.0(4) 




1.24(4) 


1.59(4) 


1.66(6) 


1.65(5) 


o- 


1.72(7) 


1.65(5) 


1.57(6) 


1.57(4) 


2++ 


0.95(5) 
0.95(2) 


1.69(6) 
1.58(4) 


1.81(7) 
2.08(15) 


1.86(11) 


2"+ 


1.87(10) 


1.83(7) 


1.79(8) 


1.76(6) 



Table 21: How the lightest SU(3) masses depend on the spatial volume. 





SU(4) ; = 28 


SU(5) ; P = 44 


state 


L = 12 


L = 16 


L = 12 


L = 16 


0++ 


1.08(2) 


1.08(2) 


1.09(3) 


1.05(2) 


++* 


1.53(6) 


1.62(4) 


1.64(6) 


1.65(3) 


o— 


1.60(8) 


1.60(4) 


1.52(7) 


1.58(5) 


2++ 


1.80(10) 


1.76(4) 


1.74(11) 


1.69(7) 


2~ + 


1.80(12) 


1.79(5) 


1.78(9) 


1.60(6) 



Table 22: How the lightest SU(4) and SU(5) masses depend on the spatial volume. 
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state 


/3 = 3.75 


/3 = 4.5 


(3 = 5 


(3 = 6 


(3 = 6 


(3 = 6 




L=8 


L=12 


L=16 


L=16 


L=24 


L=32 


++ 


2.07(9) 


1.642(43) 


-t A f~\ / - A A \ 

1.478(24) 


1.193(18) 


1.191(18) 


1.170(23) 




2.7(3) 


2.12(19) 


2.11(11) 




1.67(6) 


1.60(6) 


Q++** 






2.6(4) 




2.07(12) 


2.10(13) 


o-+ 








2.10(33) 


2.41(35) 


2.59(29) 


2++ 






2.26(17) 


1.80(8) 


1.94(11) 


1.87(11) 


2 ++* 










2.00(14) 


1.94(12) 


2-+ 






2.20(15) 


1.81(14) 


1.77(11) 


1.91(9) 


2 -+* 










2.28(27) 


2.26(16) 


1++ 








2.43(31) 


2.35(38) 


2.64(27) 


1-+ 








2.9(7) 


2.4(5) 


3.0(7) 



Table 23: The lightest SU(2) masses at lower values of (3. 



state 


P = 9 


P = 9 


P = 12 


(3 = 12 


(3 = 14.5 




L=24 


L=32 


L=32 


L=48 


L=40 


0++ 


0.7643(60) 


0.7552(67) 


0.5572(36) 


0.5628(46) 


0.4562(23) 


0++* 


1.082(16) 


1.087(14) 


0.8072(53) 


0.8047(74) 


0.6532(33) 


Q++** 


1.271(21) 


1.340(22) 


0.949(8) 


0.982(14) 


0.790(4) 


o-+ 


1.629(44) 


1.60(4) 


1.187(18) 


1.188(18) 


0.965(8) 




1.95(9) 


1.75(7) 


1.302(20) 


1.350(28) 


1.178(13) 


2++ 


1.259(13) 


1.249(15) 


0.913(7) 


0.920(13) 


0.7557(35) 


2 ++* 


1.396(22) 


1.484(27) 


0.971(8) 


1.055(11) 


0.846(7) 


2-+ 


1.276(24) 


1.286(23) 


0.928(9) 


0.910(10) 


0.7626(51) 


2 -+* 


1.532(32) 


1.480(31) 


1.089(13) 


1.094(18) 


0.958(10) 


1++ 


1.814(54) 


1.82(6) 


1.258(17) 


1.295(20) 


1.042(9) 




2.04(10) 


1.94(8) 


1.493(24) 


1.491(31) 


1.240(15) 


1-+ 


1.892(54) 


1.81(8) 


1.356(18) 


1.317(20) 


1.096(10) 


Y — I - * 


1.83(6) 


1.73(6) 


1.331(17) 


1.320(19) 


1.092(12) 



Table 24: The lightest SU(2) masses at higher values of (3. 
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state 


/? = 11 


P = 15 


P = 15 


P = 2\ 


P = 28 


P = 34 




L=12 


L=16 


L=24 


L=24 


L=32 


L=40 


0++ 


1.626(36) 


1.123(16) 


1.095(14) 


0.7561(62) 


0.5517(38) 


0.4482(36) 


o++* 


2.19(14) 


1.66(6) 


1.652(44) 


1.124(15) 


0.823(6) 


0.6737(45) 


Q++** 




2 06(9) 


2.04(11) 


1.411(23) 


1 034(9) 


8512(67) 


o — 


2 30(15) 


1 568(53) 


1 569(38) 


1.101(14) 


8133(57) 

\J • \_7 -i- '-J '-J \ l / 


6682(48) 


o — * 




2.05(16) 


2.00(10) 


1.385(21) 


1.025(10) 


0.8386(49) 


Q ** 




2.44(41) 


2.40(30) 


1.596(39) 


1.191(16) 


0.9969(85) 


(T+ 






2.32(24) 


1.627(41) 


1.206(15) 


0.9634(81) 










1 835(73) 


1 322(17) 


1.194(14) 


o+- 






2.08(23) 


1.826(59) 


1.330(15) 


1.088(10) 


0+-* 








1.98(11) 


1 582(31) 


1 315(16) 


2++ 


2.31(21) 


1 81(7) 


1.86(11) 


1.218(16) 


0.9123(70) 


0.7354(43) 


2 ++* 




2.15(14) 


2.18(15) 


1.520(27) 


1.057(10) 


0.9134(79) 


2~+ 


2.64(29) 


1.786(82) 


1.758(58) 


1.257(18) 


0.937(8) 


0.7526(48) 


2 -+* 




2.10(16) 


2.19(15) 


1.618(43) 


1.109(10) 


0.9142(65) 


2~ 




2.05(15) 


1.95(14) 


1.475(23) 


1.0928(86) 


0.8913(63) 


2 ~* 




2.21(19) 


2.35(26) 


1.705(42) 


1.254(15) 


1.0452(85) 


2+- 




2.04(16) 


1.89(11) 


1.539(27) 


1.114(11) 


0.8867(59) 


2 +-* 








1.813(54) 


1.325(15) 


1.085(11) 


1++ 








1.738(40) 


1.298(11) 


1.0513(70) 










1.933(61) 


1.481(19) 


1.212(10) 


1-+ 








1.881(51) 


1.350(12) 


1.082(8) 


j — h* 








1.902(48) 


1.352(13) 


1.096(9) 


1— 








1.780(34) 


1.269(11) 


1.036(7) 


1—* 








1.877(51) 


1.371(15) 


1.075(9) 


1+- 








1.788(37) 


1.297(11) 


1.074(9) 


j^H — * 








1.996(56) 


1.404(16) 


1.103(7) 



Table 25: The lightest SU(3) masses. 
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state 


P = 21 


P = 28 


(3 = 40 


(3 = 51 




L=12 


L=16 


L=24 


L=32 


0++ 


1.525(36) 


1.083(14) 


0.7109(52) 


0.5466(40) 


o++* 


2.31(22) 


1.616(39) 


1.080(10) 


0.821(6) 


Q++** 




1.99(8) 


1 364(21) 


1 032(9) 


o — 


2 10(13) 


1 599(36) 


1 039(13) 


8040(46) 


o — * 




2.00(7) 


1.301(18) 


1.010(8) 


Q 




2.38(17) 


1.544(24) 


1.186(10) 


o~+ 




2.35(16) 


1.575(27) 


1.200(10) 








1.80(6) 


1 325(14) 


o+- 




2.66(32) 


1.76(6) 


1.340(14) 


0+-* 






1.85(6) 


1 564(22) 


2++ 


2.08(18) 


1 76(4) 


1.168(14) 


0.9122(56) 


2 ++* 




2.09(10) 


1.408(17) 


1.085(10) 


2~+ 


2.15(19) 


1.79(5) 


1.205(18) 


0.8936(67) 


2 -+* 




1.99(10) 


1.429(23) 


1.095(12) 


2~ 




2.08(12) 


1.430(23) 


1.067(12) 


2 — * 




2.38(18) 


1.57(4) 


1.244(13) 


2+- 




2.17(10) 


1.394(24) 


1.077(9) 


2+~* 




2.53(29) 


1.59(4) 


1.296(13) 


1++ 




2.66(28) 


1.720(31) 


1.300(11) 








2.05(4) 


1.454(15) 


1-+ 




2.68(26) 


1.67(3) 


1.340(13) 


I — h* 






1.76(3) 


1.344(10) 


1 — 




2.48(13) 


1.70(3) 


1.251(9) 


I" * 






1.74(5) 


1.320(13) 


1+- 




2.42(18) 


1.73(5) 


1.264(10) 


— * 






1.79(4) 


1.313(16) 



Table 26: The lightest SU(4) masses. 
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state 


P = 33 


P = 44 


P = 64 


P = 82 




L=12 


L=16 


L=24 


L=32 


0++ 


1.550(47) 


1.054(15) 


0.695(5) 


0.5325(48) 


o++* 


2 4(3) 


1.654(30) 


1.053(11) 


0.798(7) 


Q++** 




2 06(11) 


1 342(27) 


1 009(10) 


o — 


2.29(24) 


1 581(44) 

_I_ • v_/ V7 _L \ J- J- f 


1.041(11) 


765(10) 


o — * 




2.05(11) 


1.278(20) 


0.997(8) 


Q 




2.27(19) 


1.555(34) 


1.19(2) 


o- + 




2.32(14) 


1.478(23) 


1.155(11) 








1.80(4) 


1 347(18) 


o+- 




2.54(32) 


1 78(7) 


1.256(14) 


Q+-* 






1.81(5) 


1 580(33) 


2++ 




1.69(7) 


1.08(5) 


0.8914(69) 


2 ++* 




2.05(9) 


1.39(3) 


1.072(9) 


2~+ 




1.60(6) 


1.07(5) 


0.8785(69) 


2 -+* 




2.06(14) 


1.37(3) 


1.075(9) 


2— 




2.14(11) 


1.390(25) 


1.054(10) 


2 — * 




2.57(28) 


1.559(31) 


1.225(16) 


2+- 




2.08(12) 


1.430(21) 


1.043(9 


2 +-* 




2.51(36) 


1.557(34) 


1.230(17) 


1++ 




2.38(22) 


1.718(30) 


1.294(14) 








1.93(8) 


1.421(14) 


1-+ 




2.57(23) 


1.697(31) 


1.300(10) 


J — h* 






1.78(5) 


1.331(15) 


1 — 




2.29(10) 


1.654(32) 


1.212(9) 


1 — * 






1.66(4) 


1.27(2) 






2.41(16) 


1.71(5) 


1.255(14) 


— * 






1.73(4) 


1.29(2) 



Table 27: The lightest SU(5) masses. 
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m G /y/a 




SU(2) 


STJ(3) 


STJ(4) 


SU(5) 


0++ 


4 718(43) 


4 329(41 ) 


4 236(50) 


4 184(55) 


n++* 


V . UU I X W 1 


6 52(9) 


6 38(13) 


6 20(13) 


n++** 

w 


8 15(15) 


8 23(1 7) 


8 05(22) 


7 85(22) 


o — 




6.48(9) 


6 271(95) 


6 03(18) 


* 




8.15(16) 


7.86(20) 


7.87(25) 


Q 




9.81(26) 


9.21(30) 


9.51(41) 


o-+ 


9.95(32) 


9.30(25) 


9.31(28) 


9.19(29) 


o+- 




10.52(28) 


10.35(50) 


9.43(75) 


2++ 


7.82(14) 


7.13(12) 


7.15(13) 


7.19(20) 


2 ++* 






8.51(20) 


8.59(18) 


2~+ 


7.86(14) 


7.36(11) 


6.86(18) 


7.18(16) 


2 -+* 




8.80(20) 


8.75(28) 


8.67(24) 


2~ 




8.75(17) 


8.22(32) 


8.24(21) 


2 — * 




10.31(27) 


9.91(41) 


9.79(45) 


2+- 




8.38(21) 


8.33(25) 


8.02(40) 


2 +-* 




10.51(30) 


10.64(60) 


9.97(55) 


1++ 


10.42(34) 


10.22(24) 


9.91(36) 


10.26(50) 


1-+ 


11.13(42) 


10.19(27) 


10.85(55) 


10.28(34) 


1— 




9.86(23) 


9.50(35) 


9.65(40) 


1+- 




10.41(36) 


9.70(45) 


9.93(44) 



Table 28: Glueball masses in units of the string tension: in the continuum limit. 
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best fit confidence level (%) 



state 


SU(2) 


SU(3) 


SU(4) 


SU(5) 


0++ 


85(25) 


70(20) 


70(20) 


60(20) 


0++* 


25(5) 


80(25) 


90(25) 


100(25) 


Q++** 


15(3) 


90(25) 


40(10) 


60(20) 


0~ 




95(25) 


35(10) 


15(3) 


Q — * 




90(25) 


35(10) 


23(5) 


Q ** 




65(20) 


70(20) 


70(20) 


o-+ 


90(25) 


45(15) 


85(25) 


35(10) 


o+- 




85(25) 


95(25) 


13(2) 


2++ 


95(25) 


40(10) 


25(6) 


17(3) 


2 ++* 






90(25) 


90(25) 


2~+ 


60(20) 


50(15) 


15(3) 


30(8) 


2 -+* 




12(3) 


14(10) 


65(20) 


2~ 




80(25) 


25(6) 


85(25) 


2 — * 




70(20) 


40(10) 


30(10) 


2+- 




25(5) 


45(10) 


3(0.3) 


2 +-* 




80(25) 


12(3) 


40(10) 


1++ 


60(20) 


40(10) 


100(25) 


22(4) 


1-+ 


60(20) 


60(15) 


10(2) 


100(25) 


1~ 




45(10) 


10(2) 


2(0.1) 


1+- 




25(5) 


12(3) 


15(3) 



Table 29: Confidence levels of the best fits in Table in brackets those of the fits that 
provide the errors. 
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atm with at — 0.25a, 





ft — 4 


ft — 5 3 


ft — 8 




1 2 2 60 


1 6 2 64 

J. w Ut: 


24 2 96 

in: U\J 


11 Ua lUUp 


n zm (7^ 


U .OUtJ I O J 


n i r^o\ 


0++ 


n 4.4.1(5} 


321 (A\ 


907f9'l 


D++* 
u 


U.UOI 1 J 


AAR('\ 9) 


n ^n9f^ 


n++** 
u 


fi 77(9} 

U. ( ( J 


U.OOZ I / J 


U.OUO ^ i J 


o~+ 


1.01(3) 


68f3) 


0.427(9) 


-+* 


1.14(3) 


0.73(3) 


0.527(5) 


2++ 


0.717(10) 


0.532(8) 


0.344(4) 


2++* 


0.86(2) 


0.636(9) 


0.395(6) 


2-+ 


0.750(13) 


0.538(5) 


0.339(4) 


2 -+* 


0.93(3) 


0.647(8) 


0.407(5) 


1++ 


1.03(3) 


0.75(5) 


0.471(4) 




1.3(1) 


0.76(6) 


0.551(7) 


1-+ 


1.12(4) 


0.80(6) 


0.499(6) 


I — h* 


1.22(5) 


0.80(5) 


0.477(12) 


(TrU Ps /N c ) 


0.64060(7) 


0.73195(3) 


0.82418(1) 


(TrU Pt /N c ) 


0.91436(2) 


0.93597(1) 


0.95795(1) 



Table 30: Masses with asymmetric SU(2) action; r = 0.25 in eqn([8]). Also shown are the 
average timelike and spacelike palquettes. 



p 


L s 


r 


Tl 


p = 0,l 


p = 0,2 


a t mi/a s mi 


8.0 


24 


0.25 


0.232 


0.249(11) 


0.252(5) 






16 


0.25 


0.232 


0.244(9) 


0.241(4) 


0.235(8) 


5.3 


16 


0.25 


0.221 


0.239(25) 


0.235(7) 






12 


0.25 


0.221 


0.235(10) 


0.238(4) 


0.241(6) 


4.0 


12 


0.25 


0.209 


0.236(28) 


0.241(8) 






8 


0.25 


0.209 


0.223(7) 


0.219(5) 


0.241(12) 



Table 31: Various estimates of a t /a s as described in Appendix D. Also shown is rj, the 
mean-field improved value of r. 
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state 


SU(2) ; r=1.0 


SU(2) ; r=0.25 


0++ 


4.718(43) 


4.65(10) 


0++* 


6.83(10) 


6.83(20) 


Q++** 


8.15(15) 


8.39(33) 


o-+ 


9.95(32) 


9.23(38) 


2++ 


7.82(14) 


7.81(20) 


2++* 




8.86(30) 


2-+ 


7.86(14) 


7.54(20) 


2 -+* 




8.94(27) 


1++ 


10.42(34) 


10.51(24) 


1-+ 


11.13(42) 


11.03(30) 






[10.38(44)] 


V*/g 2 


0.3353(18) 


0.3375(130) 



Table 32: Comparison between the continuum mass ratios, m G /^Ja, obtained with the 
asymmetric r=0.25 SU(2) action, and our previous r = 1 SU(2) results. Also shown is y/u/g 2 
for both cases. 
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state 


]im Nc ^oo m/g 2 N c 


slope 


N c > 


CL% 


0++ 


0.808(11) 


-0.070(79) 


2 


90(25) 


o++* 


1.227(25) 


-0.31(18) 


2 


65(20) 


Q++** 


1.581(42) 


-0.84(28) 


2 


50(15) 


o-+ 


1.787(60) 


-0.50(51) 


2 


85(25) 


2++ 


1.365(33) 


-0.25(28) 


2 


35(10) 


2-+ 


1.369(36) 


-0.20(27) 


2 


15(4) 


2 -+* 


1.704(70) 


-0.74(88) 


3 


95(25) 


1++ 


1.98(8) 


-0.90(57) 


2 


80(25) 


1-+ 


1.99(8) 


-0.61(70) 


2 


35(10) 


0~ 


1.167(42) 


0.26(50) 


3 


65(20) 


* 


1.508(72) 


-0.07(87) 


3 


65(20) 


Q 


1.77(13) 


0.24(161) 


3 


30(8) 


o+- 


1.87(23) 


0.63(245) 


3 


45(10) 


2— 


1.57(8) 


0.40(93) 


3 


55(15) 


2 ~* 


1.87(12) 


0.23(143) 


3 


90(25) 


2+- 


1.59(10) 


-0.37(117) 


3 


65(20) 


2 +-* 


1.97(17) 


-0.28(188) 


3 


50(15) 


1— 


1.85(13) 


-0.33(149) 


3 


55(15) 


1+- 


1.87(16) 


0.37(200) 


3 


45(10) 



Table 33: The large N c limit of the mass spectrum in units of g 2 N c ; with the slope of the 
linear fit when plotted against l/iV c 2 . Also the range of colours fitted and the confidence level 
of the fits. 



72 



state 


limjv c -^oo ml Jo 


slope 


N c > 


CL% 


0++ 


4 065f55) 


2 58(42) 


2 


80(25) 


Q++* 


6 18(13) 


2 68(100) 


2 


70(20) 


Q++** 


7 99(22) 


79(160) 


2 


50(15) 


0"+ 


9.02(30) 


3.52(275) 


2 


85(25) 


2++ 


6.88(16) 


3.50(134) 


2 


30(10) 


2"+ 


6.89(21) 


3.13(162) 


2 


20(5) 


2 -+* 


8.62(38) 


1.69(165) 


3 


90(25) 


1++ 


9.98(25) 


1.78(203) 


2 


80(25) 


1-+ 


10.06(40) 


3.58(365) 


2 


30(8) 


0~ 


5.91(25) 


5.24(300) 


3 


55(15) 


* 


7.63(37) 


4.61(460) 


3 


70(20) 


Q ** 


8.96(65) 


7.2(80) 


3 


35(10) 


o+- 


9.47(116) 


9.7(12.4) 


3 


40(10) 


2— 


7.89(35) 


7.6(44) 


3 


60(20) 


2 — * 


9.46(66) 


7.6(77) 


3 


95(25) 


2+- 


8.04(50) 


3.2(60) 


3 


60(20) 


2 +-* 


9.97(91) 


5.1(10.0) 


3 


50(15) 


\— 


9.36(60) 


4.4(70) 


3 


60(20) 


1+- 


9.43(75) 


8.4(98) 


3 


50(15) 



Table 34: The large N c limit of the mass spectrum in units of the string tension; with the 
slope of the linear fit when plotted against 1/iV^. Also the range of colours fitted and the 
confidence level of the fits. 
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m G /> : U(l) 


state 


(5 = 2.0 


P = 2.2 


p = 2.3 


0++ 


3.54(9) 


3.29(23) 


3.36(17) 


0~ 


1.97(7) 


1.52(5) 


1.50(5) 


o-+ 


7.1(9) 


8.1(4) 


8.8(3) 


o+- 




8.5(4) 


10.4(5) 


2++ 


5.2(7) 


4.64(30) 


5.12(27) 


2-+ 


4.8(11) 


4.9(9) 


5.45(40) 


2~ 


6.1(4) 


7.0(3) 


6.3(4) 


2+" 


7.0(5) 


6.7(3) 


6.7(6) 


1++ 




9.6(6) 


11.0(6) 


1-+ 




9.7(6) 


11.6(6) 


1— 


7.8(8) 


7.9(3) 


8.7(3) 


1+- 


8.2(5) 


8.0(3) 


8.6(3) 


a^Jo 


0.2251(18) 


0.1734(16) 


0.1505(15) 



Table 35: The U(l) mass spectrum in units of the string tension, at several values of /3. In 
the last row is the string tension in lattice units. 
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amp 




Figure 1: Mass of periodic flux loop, amp, against its length, L, at (5 — 6. The straight line 
is to guide the eye. 
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Figure 2: Mass of periodic flux loop of length, L 
correction to the linear rise. Curve is fit using eqn 



at (3 — 9; divided by L to expose the 
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Figure 3: The values of f3a^a plotted against 1//3 for SU(2). Also shown is the leading-order 
strong coupling prediction at low (3, and a leading-order continuum extrapolation at high (3. 




Figure 4: The values of (3a^a plotted against 1/(3 for SU(3). Also shown is the strong coupling 
prediction to 0{j3) at low (3, and a leading-order continuum extrapolation at high (3. 
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Figure 5: As in Fig.^j but using the mean-field improved coupling, j3j, in place of (3. 
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Figure 6: As in Fig.f| but using the mean-field improved coupling, fli, in place of f3. 




Figure 7: The value of \f& / g 2 as a function of N c . The line shows the fit in eqn fl44|) 
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Figure 8: The \ 2 P er degree of freedom against the power, a, of the leading large- N c correction 
when fitting y / a/g 2 N c . 
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Figure 9: How the lightest SU(2) masses vary with the spatial volume, L 2 , at (5 — 9. States 
are the ++ (o), the 2++(A), the 2" + (o), the l ++ (x) and the l~ + (*). Also shown is twice the 
mass of the periodic flux loop (•). 
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Figure 10: The effective masses obtained on a 24 2 96 lattice at (3 = 8 in SU(2), with a very- 
small temporal lattice spacing: a t ~ 0.25a s . States are the ++ («), the ++ *(o), the 2 ++ (*), 
the 2++*(o) and the l ++ (x). 
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Figure 11: The ratio of some SU(3) masses to ay/a, plotted against a 2 a to show how they vary 
with a: the ++ («), the (o), the 2~ + (*) and the 2 (o). Extrapolations to the continuum 
limit are shown as straight lines. 
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0.5 



0' 1 1 1 1 1 1 

0.05 0.1 0.15 0.2 0.25 

Figure 12: Some of the C = + glueball masses for 2,3,4 and 5 colours, in units of g 2 N c and 
plotted against l/iV c 2 : ++ («), ++ *(x), 2 ++ (V), 0~ + (o), l ++ (o). The best linear extrapola- 
tions to the N c = oo limit are also shown. 
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Figure 13: Some of the C = — glueball masses for 3,4 and 5 colours, in units of g 2 N c and 
plotted against l/iV^: (•), *(x), 2 (★), 1 (o). The best linear extrapolations to the 
N r = oo limit are also shown. 
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